Subject:  Acoustic  Diffraction 

Part  I-  Plane  Diffractors  and  Wedges 


Abstract:  This  memorandum  documents  the  theoretical 

investigations  in  the  Acoustic  Diffraction 
Program.  This  report  discusses  the  acoustic 
diffraction  and  backscattering  phenomena  for 
plane  and  wedge  scatterers , which  are  insonified 
by  plane  or  point  sources.  The  theories  of 
diffraction  used  in  this  report  are  those  of 
the  approximate  integral  representations  of 
Kirchof f-Rubinowicz . Those  were  compared  with 
the  geometrical  theory  of  diffraction  (GTD)  which 
is  developed  by  J.  B.  Keller,  and  is  based  on  the 
ray  theory.  It  was  found  that  the  integral 
representation  agrees  well  with  the  GTD  at  high 
frequencies  for  infinite  edges  and  long  curved 
edges  only.  The  agreement  was  worse  when  the 
scatterer  size  or  edge  length  is  finite.  New 
expressions  for  diffraction  from  finite  scat- 
terers and  edges  were  derived  and  the  deviation 
of  the  GTD  from  these  solutions  were  estimated. 
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ACOUSTIC  DIFFRACTION 

1.  INTRODUCTION 

i ^ 

’ In  dealing  with  acoustic  diffraction  it  is  expedient  to  dis- 

tinguish between  three  groups  of  problems.  For  this  reason  the 
report  will  be  divided  into  three  basic  parts  which  distinguish 
the  types  of  diffractors  : (I)  plane  diffractors  and  wedges,  (II)  three 
dimensional  bodies,  and  (III)  structures  whose  vibrations  are  coupled 
with  the  diffracted  field. 


/ 

i 


1 . 1 Plane  Diffractors  and  Wedges  - Part  I 

The  first  and  simplest  group  deals  with  plane  screens,  plane 
wedges  and  plane  apertures.  Considerable  effort  has  been  made  in 
the  past  to  deal  with  such  diffractors.  Exact  solutions  have  been 
derived  for  the  straight  edge,  the  wedge,  and  for  the  circular  and 
elliptical  apertures  and  screens,  with  both  rigid  and  pressure 
release  surfaces.  Recently,  elastic  wedges  have  also  been  in- 
vestigated and  exact  solutions  derived. 

The  exact  methods  are  very  complex,  and  while  they  yield 
numerical  solutions,  they  provide  relatively  little  insight  into 
the  physical  phenomena  causing  the  diffraction.  Simple  methods  have 
therefore  beer*  developed  ever  since  diffraction  was  first  studied. 

The  most  famous  method  is  Kirchoff's  method,  where  an  exact  surface 
integral  for  the  diffracted  field  is  solved  by  replacing  the  exact 
diffractor's  surface  boundary  conditions  by  approximate  ones.  The 
approximation  involves  the  assumptions  that  the  diffractor's  surface 
is  a totally  absorbant  one,  and  that  there  is  no  acoustic  field 
whatsoever  immediately  behind  the  diffracting  surfaces.  The  field 
in  front  of  the  screen  and  aperture  is  assumed  to  be  identical  with 
that  of  the  incident  wave.  It  is  found  by  comparison  with  exact 
solutions  that  this  heuristic  procedure  does  indeed  lead  to  very 
useful  solutions. 

Because  of  the  success  of  the  Kirchoff  theory,  solutions  have 
been  derived  by  refining  the  Kirchoff  assumptions  (W.  Braunbek,  1950, 
1954,  and  1959)  or  by  developing  a method  of  successive  approximations 
that  consists  of  pairs  of  steps,  one  satisfying  the  differential 
equation,  the  other  satisfying  the  exact  boundary  conditions 
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(W.  Franz,  1949,  1950,  1957,  Schelkunoff,  1951).  These  refined 
methods  lead  to  an  improvement  of  the  Kirchoff  solution;  however, 
since  the  deviation  from  the  exact  solution  is  still  relatively 
great  in  some  angular  range  or  the  methods  do  not  converge,  it 
hardly  seems  worthwhile  to  use  them. 

Recently,  a semi -empirical  theory  has  been  derived  by  Keller 
(1956)  using  the  concept  of  rays.  The  distance  variation  of  the 
solution  is  derived  on  the  basis  of  geometrical  optics,  on  the 
(implicit)  assumption  that  the  diffracting  edges  are  large 
compared  to  the  wave  length,  so  that  only  the  stationary  phase 
points  (see  Appendix  I)  contribute  to  the  diffraction  field.  The 
directivity  factor  then  is  deduced  by  comparison  of  the  geometrical 
optics  solutions  with  the  exact  solutions.  This  is  permissible 
since,  in  the  frequency  range  in  which  his  method  applies,  dif- 
fraction is  due  to  the  contributions  of  the  very  small  regions 
along  the  diffracting  edge  that  correspond  to  the  stationary  phase 
points  and  their  neighborhoods. 

The  geometrical  and  hence  also  the  Keller  theory  for  plane 
diffractors  breaks  down  at  caustics*  because,  in  the  geometric 
optical  limit,  the  contributions  of  a very  small  but  finite 
region  of  the  diffractor  are  condensed  at  the  caustic  to  a space 
of  zero  volume  such  as  the  space  occupied  by  a point  or  a line 
(see  pp  36) . 


* A caustic  is  a surface  at  which  adjacent  rays  intersect. 
If  this  surface  degenerates  to  a point  it  becomes  a focus.  For 
example,  for  a lens  with  spherical  surfaces  the  focal  length  de- 
pends on  the  angle  of  indicence.  The  focal  points  then  vary  with 
the  angle  of  incidence  and  if  the  incident  wave  is  spherical  so 
that  the  directions  of  the  incident  rays  vary  within  a certain 
angular  range,  the  focal  points  will  lie  on  a surface,  i.e., 
on  the  caustic  surface. 
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This  accounts  for  the  far  fetched  assumptions  needed  to 
save  this  theory.  The  geometrical  theory  also  breaks  down  whenever 
the  edges  are  of  finite  length  and  the  distance  from  the  diffractor 
exceeds  a limit  that  depends  on  the  length  of  the  edge,  the  wave 
length,  and  the  angle  of  incidence.  The  geometrical  theory  also 
breaks  down  in  practically  all  situations  when  the  diffractor  surface 
area  is  finite,  where  the  distance  from  the  diffractor  is 
greater  than  its  linear  dimensions,  and  the  angle  of  incidence 
is  not  exactly  normal. 

A theory  equivalent  to  the  geometrical  theory  of  diffraction 
for  plane  diffractors  that  are  free  of  all  the  limitations  with 
respect  to  size  of  the  diffractor  and  distance  from  it  can  be 
derived  by  exact  methods  on  the  basis  of  the  Kirchoff  assumptions. 

This  method  has  the  advantage  that  the  mathematics  is  rigorous  in 
all  its  steps,  and  that  the  derivations  are  based  on  well  defined 
assumptions.  As  a result,  we  obtain  the  correct  distance  variation 
of  the  solution  whether  the  diffractor  is  small  or  large,  and  a 
directivity  factor  that  turns  out  to  represent  a very  good  ap- 
proximation to  that  of  the  exact  solution  range,  provided  the  linear 
dimensions  of  the  diffractor  are  greater  than  a wave  length.  We 
may  then  use  Keller's  method  and  improve  the  directivity  factor 
by  comparison  with  an  exact  solution.  In  deriving  the  theory  it 
is  of  advantage  to  use  the  Kirchoff  solution  in  the  form  of  the 
Rubinowicz  integral,  (Rubinowicz,  1917)  and  to  take  the  reflected 
field  into  account  independently.  The  results  can  then  easily  be 
extended  to  surfaces  of  given  impedance. 

It  is  expected  that  the  theory  derived  in  Part  I is  adequate 
for  all  practical  computations  of  the  diffracted  field  for  plane 
diffractors,  including  the  straight  edge  and  for  sharp  wedges  of 
angular  openings  up  to  at  most  90  degrees.  However,  it  is  not  suited 
to  study  diffraction  of  blunt  bodies. 

1 . 2 Three-Dimensional  Diffractors  - Part  II 

The  diffraction  theory  for  blunt  bodies  and  for  three  dimensional 
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diffractors  in  general  is  considerably  more  complex  than  that  for 
plane  diffractors.  In  the  small  number  of  cases  that  have  been 
investigated  in  conjunction  with  this  report,  the  Kirchoff  assumptions 
seemed  to  lead  to  acceptable  results  only  when  the  linear  dimensions 
of  the  diffractor  exceeded  about  50  wavelengths.  It  has  not  yet 
been  ascertained  whether  the  Kirchoff  theory  can  be  used  as  a 
starting  point  for  improved  types  of  solutions  as  was  done  for  the 
plane  diffraction  solutions.  Results  of  such  attempts  will  be  re- 
ported separately  in  part  II. 

Keller's  method  of  deriving  the  results  represents  a considerable 
shortcut  compared  to  all  other  methods,  and  right  now  it  seems  that 
this  method  is  the  only  practical  method  to  arrive  at  useful  results 
for  the  prediction  of  the  field  produced  by  three  dimensional  dif- 
fractors. As  will  be  shown,  Keller's  expressions  can  be  obtained 
by  almost  exact  procedures.  It  will  be  shown  in  Part  II  that  the 
amplitude  variation  normal  to  the  direction  of  propagation  and  conse- 
quently also  normally  to  the  shadow  boundary  obeys  a diffraction 
equation  and  that  the  three  dimensional  Keller  theory  can  be  derived 
by  solving  this  equation  for  curved  surfaces.  Part  II  of  the 
report  will  then  represent  a basic  and  critical  review  of  the  Keller 
method  for  three  dimensional  diffractors. 

1 . 3 The  Effect  of  the  Structure  and  Structural  Vibrations  on  the 

Diffraction  - Part  III 

This  part  will  be  devoted  to  the  effect  of  the  structure  and 
of  structural  vibrations  of  the  diffractor  on  the  diffraction  field. 
Every  edge,  rib,  or  appendage  on  an  elastic  structure  represents  a 
discontinuity  of  the  vibration  field  and  leads  to  sound  radiation. 
Because  of  the  phase  correlation  of  the  vibration  with  the  incident 
wave,  the  vibration  and  diffraction  depend  on  each  other  and  cannot 
be  dealt  with  separately. 
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WHAT  IS  DIFFRACTION? 

Diffraction  is  a very  complex  phenomenon  that  can  be  interpreted 

in  a number  of  different  ways,  as  will  be  discussed  in  the  following 

sections. 

2.1  Classical  Interpretation 

The  classical  interpretation  of  diffraction  (scattering)  is 
based  on  the  ideas  of  Iluyghens.  The  diffractor  transforms  the  incident 
wave  field  into  a wave  field  that  orginiatcs  at  the  diffractor 
(scatterer)  and  that  propagates  off  into  all  directions. 

In  the  classical  theories,  scattering  and  diffraction  are 
usually  described  by  assuming  a coordinate  sys*  :m  whose  coordinate 
surfaces  are  also  surfaces  of  the  scatterer  or  diffractor.  The 
solution  is  given  by  a series  of  wave  functions  that  originate  at 
the  center  of  the  scatterer.  The  effect  of  the  scatterer  then  is 
usually  represented  by  the  acoustic  power  that  is  contained  in  the 
scattered  field.  To  obtain  this  power,  the  square  of  the  scattered 
pressure  is  computed  and  integrated  over  all  space. 

One  of  the  most  important  parameters  of  the  classical  theory 
is  the  scattering  cross  section.  It  is  defined  as  the  ratio  of  the 
scattered  power  to  the  incident  power.  Thus,  the  area  of  the  incident 
wave  front  that  transports  the  same  power  as  the  power  that  is  scat- 
tered is  called  the  scattering  cross  section.  The  classical  inter- 
pretation pretends  that  the  scattered  field  is  independent  (uncor- 
related) with  the  incident  field  and  can  thus  be  considered  as  an 
independent  entity  so  that  it  makes  sense  to  talk  about  a diffracted 
(scattered)  field  and  a diffracted  (scattered)  energy.  But  the 
fact  is  that  the  diffracted  (scattered)  field  is  not  uncorrelated 
everywhere  with  the  incident  field.  The  definition  of  the  scattering 
cross  section,  therefore,  is  somewhat  problematic. 

Incident  and  scattered  fields  are  strongly  correlated  in  the 
direction  of  propagation  of  the  incident  wave,  and  it  is  no  longer 
possible  to  consider  the  diffracted  energy  as  an  independent  entity. 

For  instance,  at  the  geometrical  shadow  boundary,  the  diffracted 
field  is  always  equal  to  half  the  incident  field  and  is  in  pliase 
I with  it  as  we  move  into  the  shadow  space  and  in  antiphase  as  we 
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move  into  the  illuminated  space.  The  scattered  field  then 
is  nothing  more  than  the  field  (the  mathematical  function)  that 
we  have  to  superimpose  on  the  incident  field  to  make  the  transi- 
tion from  the  shadow  region  to  the  illuminated  region  continuous. 
The  quantity  that  the  classical  theory  defines  as  the  scattered 
energy  is  not  related  to  any  energy  at  all,  but  represents  the 
square  of  the  deviation  of  the  total  acoustic  field  at  some 
point  from  the  incident  field,  at  that  same  point  as  if  the 
scatterer  were  absent.  However,  it  turns  out  that  as  we  move 
far  away  from  the  scatterer,  the  angular  range  near  the  geo- 
metrical shadow  boundary  over  which  incident  and  scattered 
field  are  correlated  is  small  and  decreases  to  zero  as  r 
increases  to  infinity  (for  the  wedge,  this  range  is  proportional 
to  l//r  ; see  Fig.  24.13  with  p =1,  and  0 = x/r  in  Skudrzyk, 
1971).  If  we  integrate  over  a surface  of  very  large  radius, 
the  contribution  of  the  coherence  range  can  be  neglected  and 
the  scattering  cross  section  represents  a meaningful  and  useful 
parameter  that  describes  diffraction  and  scattering. 


2.2  Diffraction  as  a Diffusion  Phenomenon 


* If  the  diffractOT  is  ideally  absorbant  , and  is  a plane 

screen,  diffraction  is  a simple  phenoirienon  that  is  basically 
connected  with  the  formation  of  a shadov/  boundary.  But  if 
the  diffractor  is  not  ideally  absorbant,  a reflected  field 

, t 

S'  is  generated.  The  reflected  field  can  be  assumed  to  be 

generated  by  an  image  source,  and  the  resultant  field  consists 
of  the  incident  waves  the  reflected  waves,  and  their  respective 
diffraction  fields.  The  exact  solutions  for  the  straight  edge 
^ ( semi -inf inite  plane)  and  for  the  wedge  are  made  up  of  two 

«.  such  fields. 
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a decomposition  of  the  solution  is  no  longer  possible  for  three 
dimensional  curved  bodies.  For  instance,  a sound  wave  that 
travels  down  a horn  is  intensified  by  multiple  reflections, 
besides  the  fact  that  the  field  that  is  generated  by  diffraction 
at  its  neck  is  also  considerably  modified  by  multiple  reflections 
and  the  diffraction  fields  connected  with  each  of  the  reflection 
fields.  Thus,  reflection  and  diffraction  can  no  longer  be 
separated . 

That  one  aspect  of  diffraction  is  due  to  the  diffusion 
of  particle  motion  into  the  shadow  spaces  follows  directly 
from  the  consideration  of  the  scattering  cross  section.  At 
very  high  frequencies,  all  the  power  that  impinges  on  the 
scattering  object  is  reflected.  One  would  thus  expect , at 
very  high  frequencies,  that  the  scattering  cross  section  is 
equal  to  the  projection  of  the  cross  section  of  the  scatterer 
on  the  wave  front  of  the  incident  wave.  But  the  computation 
shows  that,  at  very  high  frequencies,  the  scattering  cross 
section  is  twice  as  large  (J . Y.  Bowmann , T.  B.  Senior, 

P.L.  Uslenghi,  1969).  This  means  that  rays  that  have  passed 
many  wave  lengths  beyond  the  scatterer  are  still  being  scat- 
tered. Still,  in  the  geometric  optical  limit,  we  know  that 
acoustic  or  optical  coupling  takes  place  over  distances  of 
millions  of  wavelengths.  The  explanation  of  what  happens  is 
a relatively  simple  one.  The  scatterer  also  generates  a 
shadow  zone.  The  boundary  of  this  shadow  zone  is  a region 
of  discontinuity  and  consequently  acts  like  a distribution 
of  sound  sources.  These  sound  sources  emit  radiation  in  all 
directions.  At  distances  sufficiently  far  away  behind  the 
scatterer,  the  shadow  region  is  filled  in  and  the  intensity  is 
again  the  same  everywhere  as  in  the  incident  wave.  Thus  the 
incident  wave  has  to  supply  not  only  the  reflected  energy  but 
also  the  energy  that  is  needed  to  fill  in  the  shadow  region 
behind  the  scatterer.  At  very  high  frequencies,  both  phenomena 
require  the  same  amount  of  energy  and  the  scattering  cross 
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section  increases  with  the  radius  of  the  cylinder  over  which 
we  compute  the  square  of  the  diffracted  field  and  becomes 
infinite  as  this  radius  becomes  infinite.  Since  the  shadow 
space  behind  the  edge  is  infinitely  extended,  an  infinite 
amount  of  energy  is  required  to  fill  the  shadow  space  and 
consequently  the  scattering  cross  sect  ions  become  infinite. 

The  diffraction  of  energy  into  the  shadow  space  is  thus 
a consequence  of  the  diffusion  of  momentum  in  a direction 
transverse  to  the  direction  of  propagation.  The  wave  amplitudes 
vary  drastically  transverse  to  the  direction  of  propagation 
in  directions  parallel  to  the  wave  fronts.  To  satisfy  the 
wave  equation  this  variation  must  satisfy  differential  equations 
of  the  diffusion  type.  The  sound  field  diffuses  into  the  shadow 
space.  If  this  aspect  of  diffraction  is  carried  on  for  three- 
dimensional  diffractors,  we  arrive  at  solutions  similar  to 
the  Keller  solution  for  the  sound  propagation  along  curved 
surfaces . 

2 . 3 Diffraction  as  a Phenomenon  Caused  by  the  Shadow  Boundary 

The  contribution  of  the  shadow  boundary  to  the  diffraction 
phenomenon  is  of  particular  interest.  The  shadow  boundary 
represents  a surface  of  discontinuity,  - at  least  in  the  geometric 
optical  limit  - and  consequently,  according  to  Huyghens  principle 
or  in  the  Kirchoff  approximation,  acts  like  a continuous  dis- 
tribution of  dipole  sources  whose  axes  are  normal  to  the  shadow 
boundary.  For  a plane  incident  wave,  the  wave  amplitude  of  the 
incident  wave  is  constant.  In  the  geometric  optical  limit, 
the  acoustic  field  near  the  shadow  boundary  is  therefore  the 
same  everywhere. 

If  we  look  at  a homogeneous  plane  shadow  boundary  and 
apply  the  Helmholtz  Huyghens  integral  to  it  (see  Appendix  II), 
we  find  that  the  field  amplitude  decreases  exponentially  with 
the  distance  perpendicular  from  it.  In  the  case  of  a plane  incident 
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wave,  the  diffracted  wave  field  is  generated  solely  near  the 
edge  that  defines  the  shadow  boundary.  The  shadow  boundary 
itself  does  not  radiate.  If  the  incident  wave  is  cylindrical 
or  spherical,  radiation  can  also  be  generated  by  the  shadow 
boundary.  So  both,  the  region  near  the  edge,  and  the  whole 
shadow  boundary  may  radiate. 

To  investigate  the  question  whether  edge  waves  are  real, 
we  look  at  the  diffraction  of  a Dirac  pulse.  Let  the  pulse 
travel  first  as  a plane  wave  pulse,  and  let  it  be  diffracted 
at  a straight  edge.  The  computation  then  shows  (F.  Oberhettinger , 
1958)  that  this  pulsed  wave  continues  along  the  illuminated 
region  up  to  the  shadow  boundary.  The  shadow  boundary  then  does 
not  radiate  because  it  is  homogeneous  in  directions  parallel  to 
the  wave  front.  An  edge  wave  does  not  exist  except  for  points 
extremely  close  to  the  edge  or  the  shadow  boundary.  But  if 
the  incident  Dirac  pulse  has  a cylindrical  wave  front,  the 
incident  wave  amplitude  near  the  shadow  boundary  is  no  longer 
constant  but  decreases  inversely  proportionally  to  the  square 
root  of  the  distance  from  the  edge.  The  first  scattered  con- 
tribution then  arrives  at  a time  that  corresponds  to  the 
shortest  distance  from  the  source  to  the  edge  and  from  there 
to  the  field  point.  But  instead  of  a Dirac  pulse  at  only  the 
shadow  boundary,  we  observe  through  the  entire  shadow  space 
a significant  diffraction  pulse  whose  amplitude  decreases  as 

i/iT  c^t^  - const.,  i.e.  inversely  proportional  to  the  distance 
from  a point  that  travels  with  sound  velocity  in  the  direction 
of  the  incident  radiation  along  the  shadow  boundary.  The 
original  pulse  now  is  drawn  out  because  all  points  from  the 
shadow  boundary  radiate.  That  the  amplitude  decreases  not  as 

but  as  1/ct  could  have  been  expected.  The  shadow  boundary 
for  a Dirac  pulse  can  be  interpreted  as  a sequence  of  line 
elements  of  width  dx  and  of  infinite  length  (more  correctly 
as  semi  circular  rings).  But  such  a strip  of  ring  (because 
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its  radius  of  curvature  is  very  great)  acts  like  a very  long 

line  array.  Its  directivity  function  is  ^sih'^'e”  ^ ^ ^ k&>>l 

and  radiates  as  if  all  the  field  were  coming  from  a single  point 
at  its  center,  the  point  opposite  the  field  point  (E.  J.  Skudrzyk, 

1971,  page  599).  The  pressure  decreases  as  1/ct  where  ct  is 
the  distance  of  this  energized  strip  element  dx  along  the  shadow 
boundary  from  the  field  point. 

The  solution  seems  very  different  for  a train  of  sinusoidal 
oscillations.  The  combined  effect  of  edge  and  shadow  boundary 
then  becomes  equivalent  to  a wave  train  that  originates  at  the 
stationary  phase  point  (see  Appendix  I).  The  effect  of  the 
shadow  boundary  no  longer  dominates  (see  page  ) in  the  solution. 

Ke  end  up  with  Fermat's  principle  (see  Appendix  III)  which  is  more  of 
a mathematical  than  a physical  nature.  The  diffracted  field 
originates  at  and  near  a point  of  the  edge  such  that  the  sum  of 
the  distances  from  the  source  to  the  edge  point  and  from  this 
edge  point  to  the  field  point  is  a minimum  (see  page  24  and 
Appendix  IV) . The  incident  wave  field  acts  like  a continuous 
sequence  of  Dirac  pulses  of  varying  amplitudes.  For  the  first 
few  pulses  the  shadow  boundary  may  be  expected  to  play  a decisive 
role,  but  the  pulses  interfere  and  after  a very  short  time  the 
diffracted  field  seems  to  originate  from  the  stationary  phase 
points.  The  diffracted  pulse  arrives  with  a time  equal  to  the 
distance  from  the  edge,  if  the  incident  wave  is  plane,  but,  if 
the  incident  wave  is  spherical  the  pulse  arrives  on  time,  as  if 
it  had  travelled  along  the  shadow  boundary  to  a point  nearly  perpen- 
dicular to  the  field  point  and  from  there  to  the  field  point.  This 
means  that  for  sinusoidally  modulated  pulses^  the  greater  the 
curvature  of  the  wave  front  the  greater  will  be  the  distortion  of 
the  shape  of  incident  pulses  by  edge  diffraction.  These  conclusions 
can  be  easily  drawn  by  computing  the  group  velocity  delay  from  the 
solution  for  the  straight  edge. 
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3.  THE  KIRCllOFF  THEORY  OF  DIFFRACTION 

The  exact  treatment  of  diffraction  is  relatively  complicated 
because  it  must  include  multiple  reflections  as  well  as  multiple 
diffractions  and  the  reflection  of  diffracted  waves.  For  plane 
diffractors,  the  diffracted  field  decomposes  naturally  into  that 
of  the  incident  wave  and  that  of  the  reflected  wave.  But  for 
three-dimensional  diffractors,  such  a decomposition  is,  in 
general,  not  possible.  Kirchoff  simplified  the  treatment  of 
diffraction  considerably  by  assuming  all  diffractors  to  be  ideally 
abosrbent.  Primary  reflections  and  multiple  reflections  do  not 
occur  and  all  that  needs  to  be  considered  is  the  incident  wave 
and  its  diffraction  field.  This  field  can  be  easily  approximated 
and  further  improved  by  computation.  If  the  diffractor  is  plane, 
the  reflected  field  solution  can  be  determined  by  an  image  method. 
As  a consequence  of  neglecting  multiple  diffraction  and  reflection 
and  the  reflection  of  diffracted  waves,  the  Kirchoff  solution 
depends  only  on  the  shadow  boundary  and  computations  are  very 
simple . 

It  is  understandable  that  the  Kirchoff  method  will  lead  to 
good  approximations  for  plane  screens.  In  fact,  it  leads  to 
excellent  approximations.  But  the  Kirchoff  approximation  may 
fail  grossly  for  three-dimensional  diffractors.  Multiple  re- 
flection and  diffraction  can  then  no  longer  be  neglected.  In 
the  Kirchoff  solution,  a horn  shaped  diffractor  is  absorbent  on 
its  inside,  and  the  amplitude  of  the  wave  that  reaches  the  neck 
is  assumed  to  be  the  same  as  that  of  the  undisturbed  incident  wave. 
That  this  assumption  is  grossly  wrong  is  apparent. 

3 . 1 The  Exact  Huyghens -Helmholtz  Diffraction  Integral  and  the 

Kirchoff  Theory  of  Diffraction 

Probably  the  deepest  insight  into  the  theory  of  diffraction 
is  obtained  with  the  aid  of  the  (exact)  Helmholtz -Huyghens 
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integral  (see  Appendix  II)  for  the  sound  pressure  p: 


P “ T? 


IP  If 


If 


where  g is  the  infinite  free-space  Green's  function 


g = g 


(r)  - i e-f'r 


and  where  it's  derivative  in  the  n-direction  is 


-1-^- e cos  (r,n)  + a term  proportional  to  1/r^  (3) 


Integration  is  performed  over  a closed  surface  o = S + A that 
separates  the  field  point  P from  the  source  point  L,  n is  its 
outward  normal;  the  source  g(p)  = e*^^P/p  is  thus  outside  this 
surface,  and  p is  the  distance  from  the  source.  The  surfaces 
denoted  by  A and  S represent  the  aperture  and  the  screens  respect- 
ively (see  Figure  1).  In  the  above  integral,  the  first  term 
describes  the  contribution  of  the  pressure  in  terms  of  a layer 
of  dipole  sources,  because  (r,n)is  proportional  to 

cos  e * cos  (r,n).  The  second  term  is  equivalent  to  a distribution 
of  simple  sources  over  the  boundary  surfaces. 

Let  us  first  investigate  the  field  behind  the  screen  in  the 
space  T opposite  the  point  source.  If  one  represents  the  solution 
in  the  space  that  is  insonified  directly  as 

p *■  e*^^^/p  Ap  - g(p)  ♦ Ap  (4) 

where  g(p)  ■ e*^^^/p  represents  the  incident  amplitude  in  the 
light  space  and  is  zero  in  the  shadow  space  and  Ap  is  the  correction 
to  the  geometric  optical  intensity  by  the  diffraction  phenomenon. 

Kirchoff  assumes  that  p » g(p),  i.e.  that  the  pressure  is  the 
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same  as  that  in  the  incident  wave,  and  thus  neglects  any  re- 
flections and  the  effect  of  diffraction  on  the  field  values 
near  the  diffracting  object.  Kirchoff  also  assumes  p(p)  ” 0 

and  = 0 at  all  shaded  surfaces.  The  terms  containing  Ap 

will  thus  represent  the  deviation  of  the  Kirchoff  solution 
from  the  exact  solution.  Let  the  integrating  surface  now  be  a 
closed  surface  just  behind  the  screens  and  the  aperture,  so 
that  g(p)  and  ^re  zero  everywhere  on  S.  The  above  integral 

then  becomes 


47rp 


rg(p) 


+ Ap] 


do 


- [g(r)  4 g(p)  - g(pj  da 


S A 


da 


PKirch  * ■’’’^P 


C5) 


where 


Pk 


irch  ■ - h;  Juw  si  g(p)  - g(p)  ] dp  (6a) 
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is  the  Kirchoff  solution,  and  6p  is  the  deviation  of  the  Kirchoff 
solution  from  the  exact  solution: 


4tt6p  = 


S ♦ A 


'g(r) 


(Ap) 


- ^ 


da 


(6b) 


For  plane  rigid  diffractors  as  the  study  of  special  cases  shows, 
the  deviation  6p  is  usually  represented  by  fluctuations  of  the 
pressure  and  velocity  field  near  the  diffracting  edges.  For 
instance,  for  an  infinite  straight  edge,  the  distance  x between 
successive  maxima  in  the  plane  of  the  straight  edge  is  practically 
equal  to  the  sound  wave  length.  The  space  spectrum  of  these  rip- 
ples is  similar  to  that  of  a sinusoidal  pulse  of  finite  width.  It 
centers  at  the  carrier  wave  length  which  corresponds  to  the 
coincidence  wave  length  for  sound  radiation.  The  sound  components 
whose  wavelength  is  slightly  longer,  radiate  in  a direction  at 
or  close  to  grazing  to  the  plane  of  the  aperture  and  thus  affect 
the  sound  radiation  in  directions  normal  to  the  shadow  boundary 
whereas  the  long  wave  number  components  radiate  in  the  forward 
and  backward  direction  (see  Theory  of  Sound  Vibration  from 
Infinite  Plate,  for  instance,  E.  J.  Skudrzyk,  1971,  page  319). 

The  Kirchoff  diffraction  solution^which  assumes  ideally 
absorbent  screens  and  consequently  neglects  the  diffraction 
due  to  the  reflected  field,  vanishes  for  observation  points  which 
are  located  in  the  direction  opposite  to  the  direction  of 
propagation  of  the  incident  rays  (as  for  an  infinite  plane  wave, 
see  page  46  ).  This  conclusion  follows  from  the  Rubinowicz 

ingegral , page  22  (because  cos(p,n)  = 0).  Since  there  is  no 
Kirchoff  error^ 6p  is  responsible  for  all  the  radiation  in  this 
direction . 

If  the  screen  or  diffractor  is  of  the  pressure  release  type, 
an  additional  source  field  will  be  generated  near  the  edges  of  the 
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diffractor,  because  of  periodic  compression  of  the  material 
by  the  sound  vavo  . The  pressure  scattered  by  a small  particle 
is  given  by  RaylcTgl  , Theory  of  Sound,  Volume  II  or  Skudrzyk  1971, 
pages  416,  504*. 


r 


T 


cos 


■ 

9 


(7) 


T = volume  of  scatter 
= bulk  modulus 


p = density 

where  0 is  the  angle  between  the  direction  of  propagation  and 
the  vector  from  the  scatterer  to  the  field  (see  Figure  2)  . 

Letters  with  no  subscripts  refer  to  the  undistribued  medium, 
letters  with  subscripts  zero  refer  to  the  scatter.  Backscatter ing 
occurs  for  9 = ir . If  X^^  “ ~ i.e.  pressure  release,  we  have 


T [1  + cos  9] 


(8) 


A similar  expression  is  obtained  for  a rigid  scatterer  except 
that  the  expression  in  the  bracket  is  replaced  by  [-1  + cos  0]. 
Thus  fur  a pressure  release  scatterer  there  is  practically  no 
scattering  in  the  direction  back  to  the  source,  and  considerable 
scattering  in  the  forward  direction  where  exactly  the  opposite 
is  true  for  a pressure  release  scatterer.  We  can  expect  that 
the  Kirchoff  solution  for  scattering  along  the  axis  of  a pressure 
release  screen  (at  low  frequencies)  deviates  more  from  that  of 
the  exact  solution  than  it  does  for  a rigid  screen.  Both  the 
rigid  and  pressure  release  Kirchoff  solutions  will  be  less 
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accurate  for  scattered  rays  that  are  incident  at  or  near  grazing 
to  the  screen. 

Figure  3 shows  a comparison  between  the  exact  solution  and 
the  Kirchoff  solution  for  the  diffraction  of  a rigid  disk  for 
normally  incident  sound.  Agreement  is  very  good  for  large 
values  of  ka  up  to  angles  close  to  the  first  zero  in  the  Kirchoff 
solution.  In  the  exact  solution  the  zeroes  of  the  Kirchoff 
solution  are  filled  in,  and  scattering  for  grazing  incident  rays 
is  greater.  However,  over  quite  a wide  angular  range,  i.e. 
over  the  range  where  scattering  or  diffraction  is  significant, 
the  Kirchoff  solution  gives  a very  good  description  of  the 
phenomenon,  see  Leitner  1949. 

We  may  expect  that  at  high  frequencies,  the  velocity  over 
the  aperture  is  very  nearly  that  in  the  incident  wave  and  de- 
viations from  it  will  only  occur  near  the  edge  of  the  aperture. 

If  also  the  ratio  of  the  radius  of  curvature  to  the  wave  length 
is  great,  then  the  field  in  the  proximity  of  the  edge  will 
be  similar  to  that  near  the  edge  of  a straight  edge.  We  can 
then  use  this  field  to  correct  the  Kirchoff  solution.  W.  Braunbeck, 
in  a series  of  papers  from  1950  onwards,  developed  this  method 
using  the  exact  Sommerfeld  solution  for  the  correction  field.  It 
happens  that  this  relatively  complex  solution  when  substituted 
in  the  Kirchoff  formula  leads  to  a known  integral.  The  Braunbeck- 
Kirchoff  solution  then  leads  to  a solution  for  the  disc  or  aperture 
that  agrees  fairly  closely  with  the  exact  results  of  Meixner  1949  even 
at  relatively  low  frequencies  as  is  illustrated  by  Figure  4. 

We  thus  arrive  at  the  conclusion  that  the  Kirchoff  theory 
is  a good  approximation  for  the  phenomena  of  diffraction  by  plane 
diffractors.  It  has  the  great  advantage  of  being  based  on 
assumptions  that  are  clearly  specified  and  is  exact  in  the 
derivations  that  are  derived  from  these  assumptions. 

The  Kirchoff  theory  does  not  include  the  field  that  is 
generated  by  the  reflected  wave  nor  by  multiply  reflected  waves. 

For  plane  diffractors,  the  field  of  the  reflected  wave  can 
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easily  be  supplied  as  the  field  generated  by  the  image  source. 

But  multiple  reflections  as  it  occurs  with  three-dimensional 
diffractors  is  not  in  the  reach  of  the  Kirchoff  theory  (see  the 
case  of  the  open  horn  discussed  above),  because  all  the  dif- 
fracting surfaces  are  assumed  to  be  ideally  absorbent.  For 
the  reflected,  as  well  as  for  the  diffracted  waves,  the  dif- 
fractor  is  ideally  transparent  and  multiply  reflected  waves  are 
excluded  by  the  definition  of  the  Kirchoff  diffractor.  From 
the  point  of  the  diffracted  wave,  the  diffractor  seems  to  have 
no  volume  and  the  diffracted  waves  always  propagate  into  the 
full  angular  space  4tt.  In  contrast,  the  exact  solution  describes 
a diffracted  field  that  is  confined  to  the  space  outside  the 
diffractor.  Thus  for  a blunt  wedge,  the  free  space  around  the 
edge  extends  over  an  angle  of  only  a little  more  than  tt  radians, 
whereas  for  the  Kirchoff  space  this  angle  is  always  2n  radians. 

A ray  that  has  been  diffracted  once  is  immune  to  further  inter- 
action with  the  diffractor.  It  traverses  it  as  if  the  diffractor 
is  transparent.  Because  of  this  assumption,  the  shape  of  the 
diffractor  has  no  effect  on  the  field  near  to  it.  Therefore,  we  can 
not  expect  the  Kirchoff  tlieory  to  yield  approximations  for  the 
diffraction  field  produced  by  three-dimensional  diffractors. 

3 . 2 The  Kirchoff  Solution 

The  Kirchoff  Solution  leads  to  excellent  agreement  with  the 
exact  solution  whenever  the  direction  of  the  diffracted  rays  do 
not  differ  greatly  from  that  of  the  incident  beam,  and  it  is  a 
relatively  easy  task  to  derive  this  solution.  It  leads  to 
fair  agreement  also  for  other  angles  of  diffraction.  It  gives 
insight  into  the  irfiysics  of  diffraction  like  no  other  method  or 
type  of  solution  does. 

Because  the  Kirchoff  solution  agrees  closely  with  the 
exact  solution  in  certain  angular  ranges,  it  is  a theory  that 
is  well  suited  to  derive  the  distance  dependence  of  the  diffracted 
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field.  It  is  very  suited  as  nuclus  for  an  improved  solution,  for 

instance,  by  applying  the  Keller  method  to  the  Kirchoff  solution 

as  a starting  point  rather  than  to  the  geometric  optical  divergence 

relations.  This  procedure  will  safeguard  us  from  comitting  any 

serious  error  because  the  results  are  exact  conclusions  derived 

from  well  defined  and  very  reasonable  assumptions. 

The  Helmholtz  Huyghens  integral  (Appendix  II)  expresses  the 

field  at  a point  at  an  origin  P (r  = 0)  inside  a closed  surface 

9 d) 

in  terms  of  the  values  of  (b  and  ^ at  that  surface: 


o 


where  n is  the  outward  normal  to  the  surface  a and 

g(r)  * e’^^’"/r  (10) 


is  the  full  space  Green's  function  and  r is  the  distance  from 
the  field  point  to  the  surface.  The  sources  arc  outside  the 
surface  a,  so  that  o separates  the  sources  from  the  field  point. 

If  the  field  point  and  the  sources  are  on  the  same  side  of  the 
surface  (either  outside  or  inside)  then  the  above  integral  is  zero 


(11) 


where 


g ”=  e'-’'*^^/r  (12) 

and  r is  the  distance  from  any  fixed  field  point  on  the  source 
side  to  the  surface  of  integration. 

It  is  convenient  to  separate  the  geometric  optical  field 
from  the  diffraction  field.  This  is  done  by  considering  the 


i 
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volume  outside  the  volume  taken  up  by  the  screen  as  integration 
volume  (see  Fig,  5)  and  both  sides  of  the  screens  as  tlie  surface 
of  integration.  The  source  now  is  contained  in  the  space  of 
integration  and  must  be  excluded  by  a little  sphere  in  the 
same  manner,  as  was  excluded  the  immediate  vicinity  of  the 
field  point  previously  from  the  space  of  integration.  Source 
and  field  points  thus  are  surrounded  by  little  spheres  and 
(except  for  the  sign)  contribute  in  the  same  manner. 

The  integral  over  the  source  sphere  then  is  equal  to  the 
volume  flow  of  the  source  (unity  here)  taken  with  a sign 
opposite  to  that  of  the  surface  integral  around  the  field  point 

I 

and  the  source  term  thus  appears  on  the  right  hand  side  of 

with  positive  sign,  n is  the  normal  pointing  into  the  screen. 

If  the  screens  are  thin,  we  can  combine  the  contribution  of  I 

I the  illuminated  and  the  shadow  side  of  the  screens  and  integrate 

: only  over  the  surface  of  the  screen. 

If  we  do  this  with  the  Helmholtz  Huyghens  integral,  Equat ion  (8), 
we  obtain  the  result 

screen 
ilium,  side 

The  normal  now  points  into  the  illuminated  space.  Here  we 

cJ  6 

have  an  exact  result,  and  4>  and  ^ represent  the  discontinuities 

in  the  pressure  and  pressure  gradient  (velocity)  across  the 
screen  surfaces. 

The  Kirchoff  integral  is  obtained  by  replacing  i)  and  ^ 

by  their  values  in  the  incident  wave  as  if  the  screens  were 
^ perfectly  abosrbcnt,  and  assuming  0 = 0 , = 0 on  the  sliadow 

side  of  the  screens. 

The  last  assumption  means  that  0 is  identically  zero  everywhere; 

1 

I 
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this  amounts  to  neglecting  the  diffracted  field  adjacent  to  the 
screen  in  the  boundary  conditions  for  the  shadow  region. 

The  incident  wave  can  be  represented  by  an  exactly  similar 
integral  over  a closed  surface  (screen  + aperture) , The  normal 
n then  is  positive  from  the  shadow  space  to  the  light  space 
(towards  the  source).  In  the  above  integral,  n'  is  chosen 
positive  in  the  opposite  direction  (n'  = - n)  . By  changing  ii' 
into  -it,  and  replacing  the  incident  wave  by  the  surface  integral, 
the  parts  over  the  screen  have  opposite  signs  and  we  are  left 
with 

* ■ |g(r)  - g(p)  g(r)j  da 

ilium,  screen  ' 

+ aperture 

|g(r)  . g (p)  i|Ill  I da 

ilium,  screen  ■' 


+ 


1 


9g(p)  . 

8n 


g(p) 


an 


g(r) 


da 


aperture 


(14) 


The  last  formula  represents  the  deviation  of  the  Kirchoff  formula 
from  the  exact  solution. 

Next  we  introduce  the  integral  over  the  aperture  A and  the 

truncated  shadow  (i.e.  the  shadow)  K,  (see  Fig.  6).  This  integral 

- i kR 

obviously  represents  the  incident  wave  g(R)  = e /4ttR 


♦g  = g(R) 


represents  the  illumination  of  the  aperture  if  diffraction  is 
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neglected.  The  integral  we  are  interested  in  (Eq.  14)  is  exactly 
as  the  one  above  (Eq.  15)  except  that  it  is  taken  only  over  the 
surface  A.  Thus 


4> 


K 


$ . $ 
g cone 


of 


This  integral  represents 
the  incident  field  and 


the  resultant  field, 
the  diffraction  field 


i . e . the 
Hence 


sum 


K g 


(j  = (j)  + 

cone  g 


(16) 


The  last  form  of  the  right  band  side  defines  the  diffracted 
field 


4>j  = - 

d cone 


(17) 


We  have  thus  transformed  the  Kirchoff  integral  into  an  integral 
over  the  shadow  boundary  K (K  = shadow  cone  surface:) 


g(r) 


3 


g(p)  - g(p)  g(r)>  da 


(18) 


This  integral  now  describes  the  diffraction  as  the  exclusive  effect 
of  the  shadow  boundary.  In  the  Kirchoff  approximation  the  J:.f- 
fractors  that  cast  the  same  geometrical  shadow  generate  the  same 
diffraction  field. 


4.  THE  RUBIN'OWICZ  IN’TEGR..\L 

The  Kirchoff  integral  (Eq . 18)  over  the  shadow  boundary  can 
be  transformed  with  a line  integral  over  the  diffracting  edge. 
The  surface  p = constant  (p  is  the  distance  along  light  ray  from 
source  point)  intersects  with  the  surface  K at  right  angles. 
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Hence  at  K we  must  have: 


V 


n 


d 

Fir 


g(p) 


= 0 
K 


(19) 


Therefore,  the  diffracted  field  becomes 


' W ij  g(p)  g(r)  do  = 

K K 


jk(p  + r) 
pr 


(1  + (Ti  ,r)do 


(20) 


where, as  before, p is  the  distance  from  the  source  and  r that  from  the 
field  point  to  the  element  do  of  the  surface  of  integration.  Thus, 
every  element  of  the  shadow  boundary  acts  like  a dipole  source, 
whose  amplitude  is  proportional  to  that  of  the  incident  wave, 
and  whose  axis  is  normal  to  the  shadow  boundary,  where  the 
normal  vector  points  into  the  shadow  region  . It  is  possible  to 
transform  this  surface  integral  into  a line  integral  over  the 
diffracting  edge  by  performing  a somewhat  tricky  partial  integra- 
tion (see  Appendix  IV) . The  result  is 


at. 


sin  (p,ds)  ds; 


edge 


cos  (r,n) 

1 cos  (r,p) 


(21) 


L: 

< 

■ t 

H 

i 

I 


Note  that  r and  p are  vectors  that  point  to  the  edge  element  ds  , and 
that  ^ T,p  is  the  angle  subtended  between  the  positive  direction 
of  these  end  vectors  (see  Fig. 7 ).  The  factor  a is  equal  to 

unity  if  the  incident  wave  reaches  the  field  point  unobstructed 
and  zero  otherwise.  The  terms  in  the  bracket  describe  the  field 
that  is  incident  on  the  element  ds , sin  (p  ,2s)  being  the  length 
of  the  projection  of  ds  onto  the  incident  ray.  The  second  group 
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of  factors  represents  the  elementary,  but  ficticious  wave  that 
is  emitted  by  the  edge  element  ds . This  ficticious  wave  has  the 
directivity  factor  cos  0 = cos  (r,n)  as  if  it  came  from  a dipole 
with  the  axis  normal  to  the  shadow  boundary;  but  its  amplitude 
depends  in  a complex  manner  on  the  distance.  As  far  as  its 
amplitude  is  concerned,  it  seems  to  come  from  a point  at  a 
distance  r'  from  the  field  point  where 

-V  -> 

r ' = r [ 1 + cos  (r  , p)  ] , 


It  is  not  possible  to  interpret  the  whole  function  cos  (n,r)/ 
[1  + cos  (r  ,p) J as  adirectivity  factor  of  a wave,  since  this 
function  becomes  infinite  at  the  shadow  boundary  (■^r,p  = it). 
If  interpreted  as  a directivity  function,  it  would  have  to  be 
represented  by  a kind  of  delta  function.  For  a plane  incident 
wave,  and 


♦ 


1 

Ti" 


sin  (J,3s)  ds  (23) 

1 + COS  Cl'.p) 


edge 

where  the  x axis  has  been  assumed  parallel  to  the  direction  of 
the  incident  rays,  and  x is  the  x coordinate  of  the  edge  element 
ds  . 


Every  point  of  the  diffracting  edge  can  then  be  identified 
with  a sound  source  and  the  resultant  diffraction  field  is  given 
by  an  integral  of  the  form 


^d 


edge 


A(pg ,3s) 


ds 


(24) 


where  p^  is  the  distance  from  the  source  from  the  diffracting 
element,  r^  that  of  the  field  point  to  the  edge  element  ds  and 
A (Pg,r^,3s)  ds  represents  the  amplitude  of  the  edge  wave  that 
is  generated  at  the  element  ds . Unless  the  edge  has  sharp  corners, 
the  function  A will  be  a slowly  varying  function  of  s,  and  if  we 
assume  that  the  wavelength  2TT/k  is  sufficiently  large  compared  to 
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the  length  of  the  edge,  the  integral  can  be  evaluated  by  the 
stationary  phase  method.  The  stationary  phase  point  is  given 
by  the  zero  of  the  derivative  of  the  exponent  (see  Appendix  I) 

jf-lKp,  - r^]  - 0 


or 


3?  (Ps  • rp  = 0 . (25) 

This  means  that  the  path  of  the  diffracted  ray  from  the  source 
to  the  field  point  is  stationary  (i.e.  is  a maximum,  a minimum, 
or  represents  a point  of  inflection  when  plotted  as  a function 
of  s . Relation  (25)  is  identical  with  Fermat's  principle  which 
is  derived  by  assuming  geometric  optical  conditions. 

The  law  of  edge  refraction  shows  that  Fermat's  geometric 
optical  principle  applies  to  the  stationary  phase  points  of  the 
diffraction  integral  whenever  it  is  possible  to  decompose  the 
resultant  field  into  a geometric  optical  field  and  a diffraction 
field,  i.e.  whenever  Huyghens  principle,  the  Kirchoff  theory, 
or  the  Rayleigh  diffraction  integral,  or  a similar  approximation 
applies . 

The  edge  law  can  also  be  written  in  the  following  form 
d r dp 

— jg-  + — = 0 = cosCf^,^s)  + cos  (Pg,3s)  , (26) 

which  implies  that  the  stationary  phase  contribution  arises  from 
the  geometrical  relation 

cos(rg,'Ss)  = - cos  (Pg,^s)  (27) 

The  condition  that  radiation  is  emitted  to  the  field  point  is  that 
^ (p'^g.cfs)  - 7T  - ^ (r*g,d^) 


(28) 
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The  incident  ray  hits  the  edge  under  the  same  (grazing)  angle 
as  the  diffracted  ray  leaves  it  (see  Figure  8).  The  basic 
diffraction  law  is  thus  very  similar  to  the  reflection  law  of 
a mirror 


Angle  of  incidence  (with  respect  to  the  edge)  = Angle  of  Reflection 


The  only  difference  is  that  because  of  the  one-dimensionality  of 
the  edge  diffractor  the  incident  and  reflected  rays  lie  not  in 
a single  plane  as  they  do  for  a flat  mirror;  the  reflected  ray 
lies  on  the  surface  of  a cone  whose  semi-apex  angle  is  equal  to 
the  angle  of  reflection  (Fig.  8r  ).  We  shall  show  below  that 
the  amplitude  of  the  diffracted  field  is  a maximum  in  the 
direction  of  the  incident  ray,  i.e.,  at  the  shadow  boundary,  for 
which  it  is  equal  to  one  half  the  incident  amplitude,  and  that 
it  is  approximately  inversely  proportional  to  the  angular  distance 
from  the  shadow  boundary  (Fig.  "Bb)  . This  means  that  the  diffracted 
field  concentrates  on  the  surface  of  a cone,  whose  angle  at  the 
apex  with  ds  is  the  same  as  that  of  the  incident  ray.  Thus  incident 
and  diffracted  rays  lie  at  opposite  sides  of  the  diffracting  edge 
element  ds . If  we  are  interested  in  back  scattering  only,  the 
diffracted  and  incident  ray  must  be  coincident.  For  backscattering 
8=  90°,  i.e.  the  radiation  must  impinge  normally  to  the  scattering 
edge . 

The  locations  of  the  stationary  phase  points  can  easily  be 
determined  by  the  edge  law  of  diffraction , (Eq . 25).  Their 
contribution  to  the  diffracted  field  then  is  given  by  the  standard 
formula 


edge 


S (s)  e 


ug(s) 


ds 


■I 


S(s  ) /T? 


fu  g’’(sj) 


1/2 


(29) 


and  in  terms  of  the  variables  of  diffraction,  writing  ug(s)  = 
-jkc  * *jk(p  + r)  , we  have  for  each  stationary  phase  point  the 


26 


May  14,  1973 
EJS:SlH:ADS;lms 


following  contribution  to  the  diffracted  field 

^ — i—  — e*^’ 1 cos(n,r  ) ^ 

V V V V 

where  ds^  represents  an  element  of  the  edge  for  which 
^ ^ (Ps  * ’’s^  " ° • 


The  derivative  of  r^  and  with  respect  to  s are  computed  with 
the  aid  of  the  well  known  Frenet  formulaa.  We  find  (see  Appendix  V ) 


■V  = 


9*- 

9s' 


Crj-Pj) 


Sin  (p 


,as,(^.  1 


^ [cos  (Pg,h)+cos  (r^  ,h)  ] 
K 


If  this  is  substituted  above,  we  obtain  the  exact  Kirchoff  con- 
tribution of  a stationary  phase  point  on  an  edge  of  arbitrary 
radius  of  curvature  Rj,  , see  Figure  (9)  : 


1 1 
2/2i:k  p 


,-jk(p  + r)±jTT/4 


r _ / 1 . 1 . 1 fcos(p.h)  + cos(r,h)l  11  + 

Vp  f 

s ^ 


cos  I'n  ,r  1 
cos (r  ,p) 


(32) 


Where  r and  p denote  the  distances  to  the  stationary  phase  point  s^ 

as  determined  by  Eq.  25,  and  li  is  the  unit  vector  along  the  main 
normal  of  the  edge  curve,  and  is  the  radius  of  curvature  of  the 
edge.  R|^  is  basically  positive  and  points  towards  the  center 
of  curvature.  .Note  that  r points  from  the  field  point  and  p from 
the  source  point  to  the  edge  (see  Fig.  7).  The  effect  of  the 
curvature  Rj^  of  the  edge  curve  is  interesting,  as  can  be  exhibited 
from  the  following  example:  Consider  a disk  of  radius  a,  such  that 

a polar  coord  system  (r,9,4»)  is  centered  at  its  center  with  its 
axis  normal  to  the  plane  of  the  disc.  Let  the  source  be  located 
at  (i,6,o)  and  the  receiver  be  located  at  (Jl  '^)  . Thus,  the 
diffracted  field  has  two  stationary  phase  points,  s,  at  (a,T-,o) 
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P 


i : 


(minimum)  and 
first  point 
as  follows  if 


S2  at  (a,  j • tt)  (maximum)  . The 
to  the  diffracted  field  can  be 
l>>a: 


contribution  of  the 
computed  from  (32) 


^ ^I’^l  " I ^ ’^I'^l  ° ^ Pj.ds  = j 

^ n,r^  = TT  -9,  ^ ^i>Pi  “ 7 ® » r = £ + a,  p = £-a,  Rj.  = a 


Thus 


(Vi 


1 

2/TnYZ 


g-jk[2£  + a(l-sin0)]+j7T/4 
£./  ^(1  -s inG) 


cosG 

l-sinQ 


The  contribution  of  the  second  stationary  phase  = a,  tt 
is  the  same  except  that  K and  n must  be  replaced  by  -K  and  -n 
and  the  angles  relative  to  h and  n by  angles  that  are  measured 
by  ITT . 

The  contribution  of  the  second  stationary  phase  point  is 
then  given  by 


1 


I 


-jk[2£-a(l-sin0)]  -j7r/4 


2/nrn: 


l/~^  (l-sin0) 


COS0 
1 -s in0 


The  factor  l//i/a  = describes  the  strong  divergence  of 

the  diffracted  field  due  to  the  curvature  of  the  edge.  In 
performing  the  computation,  we  have  implicitly  assumed  the  existance 
of  a shadow  boundary.  For  6 = a shadow  boundary  would  only  be 

formed  if  the  disc  was  infinately  large.  If  9 = ^ and  the  disc 

was  finite,  waves  would  pass  over  both  its  surfaces  undisturbed, 
and  consequently  no  diffraction  would  be  observed. 

The  stationary  phase  evaluation  makes  sense  only  if  the 
diffracting  edge  is  sufficiently  long  relative  to  the  wave  length 
so  that  the  limits  of  integration  along  the  edge  can  be  replaced 
by  minus  and  plus  infinity.  This  requirement  is  usually  not 
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fulfilled  for  finite  edges  if  the  incident  wave  is  oblique  or  if 
I the  edges  are  short.  Thus,  a stationary  phase  point  is  either 

beyond  the  edge  and  consequently  integration  does  not  extend 
to  this  point,  or  the  integrand  may  vary  so  little  during  the 
integration  that  the  effects  of  phase  cancellation  are  neglibible, 
and  that  every  element  ds  along  the  whole  range  contributes 
about  the  same  amount  to  the  resultant  field.  The  stationary 
phase  method  leads  to  the  same  variation  of  the  solution  with 
the  distance  from  the  diffractor  of  Formats'  principle  or  the 
geometrical  theory  of  diffraction.  This  result  shows  that  the 
geometrical  theory  of  diffraction  is  applicable  to  real  finite 
diffractors  only  in  very  special  instances.  More  about  the 
solution  in  the  distance,  angle  and  frequency  range  where  the 
geometrical  theory  fails  will  follow  in  sections  10,  11  and  12. 

S . THE  GEOMETRICAL  THEORY  OF  PI FFR.\CTION 

The  geometrical  theory  of  diffraction  as  it  is  understood 
today  is  a theory  that  is  based  on  the  laws  of  geometrical  optics. 
The  missing  details  are  impressed  on  this  theory  by  comparison 
with  exact  solutions.  The  diffraction  field  obtained  from  the 
geometrical  theory  of  diffraction  has  the  same  amplitude  distance 
dependence  as  does  the  high  frequency  limit  of  the  Kirchoff- 
Rubinowicz  theory  when  diffraction  reduces  to  the  energy  emitted 
by  a single  or  a finite  number  of  stationary  phase  points. 

However,  the  geometric  optical  theory  gives  no  indication  of  the 
amplitude  of  the  diffracted  field  and  its  variation  with  angle. 
Keller  supplied  the  two  missing  items  by  replacing  his  constants 
by  angle  functionsthat  make  the  geometrical  solution  agree  with 
the  exact  solution.  In  the  high  frequency  limit,  all  plane 
diffractors  are  equivalent  to  semi - inf inite  planes.  Furthermore, 
since  the  diffraction  fields  are  generated  by  only  the  edge  elements 
near  and  at  the  stationary  phase  points,  and  since  the  divergence 
of  the  diffracted  field  is  only  a function  of  the  curvature  of 
the  edge,  it  is  thus  understandable  why  Keller's  method  of 
geometrical  diffraction  is  very  successful. 

j 
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Let  us  now  turn  to  the  derivation  of  the  basic  relations 

that  lead  to  the  geometriral  theory  of  diffraction. 

Let  us  first  consider  the  wave  front  of  a bundle  of  diffracted 

rays.  The  rays  are  the  normal  trajectories  of  the  wave  front. 

Let  us  then  consider  a tube  of  diverging  rays,  whose  corss- 

scction  is  rectangular,  with  dimensions  ds^^  and  ds2  and 

consider  a wave  front  in  the  plane  in  which  its  radii  of  curvature 

are  R and  R, . The  height  and  breadth  of  the  corss  section  of 
a 3 

the  tube  (Fig. 10)  then  are  ds , = R d9,  and  ds , = Rq’.Ot 

lal  262 


dOj  = dSj^ds2 


\^6  d0j  d02 


Let  the  cross-sectional  surfaces  of  the  tube  of  rays  represent 
surfaces  of  principle  curvatures.  Consider  an  edge  element  ds 

and  let  R^  be  a principle  radius  of  curvature  of  the  wave  front 
in  the  plane  of  dSj^  and  let  dsj  be  the  height  of  the  cross  section 


at  R * Z.  We  then  have 
a 


a?Y  = 


(^  + R ) de, 

^ OL^  1 


(£  + R ) 
^ a 


1 do, 

a 1 


A similar  relation  applies  for  the  width  of  the  tube 
ds;  ^ Rfl 


The  change  in  cross-sectional  area  of  the  tube  over  the  distance 
I is  thus  given  by 


ds|ds^ 


l"‘^2 


iZ  * R^)(Jl  + Rg) 


We  assume  that  the  power  dX  transmitted  through  an  area  da  by  a 
wave  incident  in  the  direction  of  its  normal  is  proportional 
to  the  area  of  da  and  to  the  square  of  the  amplitude  A of  the 
wave.  Since  the  same  power  passes  through  dOj  and  do2,  then 

dN  - dOjAj  - d02A2, 
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and  the  wave  amplitudes  are  related  as  follows 
, ^ da 

1 _ 

■R 


'1 


35-2 


A? 


(38) 


If  (ds)  is  the  elementary  length  of  the  scatterer,  waves  will 
propagate  outward  from  (ds)  in  all  directions.  This  means  that 
one  principal  radius  of  curvature  will  be  the  distance  from  the 
edge  (ds)  regardless  of  the  second  principal  radius  of  curvature 
of  the  incident  wave  front.  Hence  we  set  Ro<<£. 

p 


The  preceeding  equation  then  reduces  to 


A 


2 

2 


TF^Tiyi— 


The  quantity 


2nRgAj  2 


2TTconst . 


(39) 


(40) 


is  the  energy  tliat  is  emitted  per  unit  length  of  the  edge, 

2 

where  is  a constant  that  represents  the  energy  that  crosses 

the  surface  of  a cylinder  of  a very  small  radius  Rg  per  unit 
length.  Introducing  into  Eq . (39)  then  leads  to 


R 

a 

TJ; 


ITT 


and 


A-> 


o V TT-+~T7TT1T 


(41) 


If  A2  is  an  amplitude,  B^  is  an  amplitude  multiplied  by  the 
square  root  of  a length.  The  radius  of  curvature  R^  of  an 

elementary  area  do  normal  to  the  diffracted  rays  (i.e.  the 
divergence  of  the  elementary  rays)  is  determined  by  the  inter- 
section of  the  diffraction  cones  that  correspond  to  two  points 
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on  the  edge,  a small  distance  (ds)  apart  (see  Figs.  ICand  11). 
The  scattering  cones  intersect.  Let  x(s)  be  the  vector  from 
an  arbitrary  origin  to  the  scattering  element  ds  where  (s) 
is  the  distance  along  the  curved  edge.  First,  let  us  derive 
the  equation  of  a single  scattering  cone  with  the  apex  on  the 
edge  at  x(s)  and  the  half  angle  6(s)  = ^ (p,ds)  where  p is 
vector  from  the  source  to  the  edge.  The  direction  of  the  cone 
axis  is  the  same  as  that  of  the  edge  element  ds. 

A variable  point  y on  a cone  surface  v/hose  apex  is  at 
x(s)  satisfies  the  equation 


[y  - x(s)]  » x(s) 
|y  - x(s) I 


COS0  (s)  , x(s) 


3x(s)  ds 

“Ti — a?  ’ 


i(s)|  = 1. 


(42) 


A second  cone  is  obtained  by  increasing  s by  ds . The  curve 
of  intersection  of  two  cones  with  the  apexes  at  x(s)  and 
x(s  + ds)  is  given  by  eliminating  s from  the  two  equations. 
This  is  best  done  by  retaining  Fq . (42)  and, instead  of  writing 
out  the  equation  with  the  apex  at  x(s  + ds),and  forming  the 
difference  between  the  two  equations,  one  can  differentiate 
Eq.  42  with  respect  so  s, 

[y  - x(s)  ] »x(s) -x(s)  *x(s)  =- |y-x(s)|  sinB(s) 


. *x(s)  cos6(s) 

|y  - x(s)  I 


(43) 


Because  of  Eq . (42),  the  last  term  on  the  right  reduces  to 
-cos  6(s) . If  the  square  root  in  the  first  term  on  the  right 

is  represented  by  Eq.  (42)  and  with  the  aid  of  the  relations: 

• • • • 

x*x  ■ |x(^  ” 1,  x(s)  * 5i/Rj^,  x(s)  t(s),  we  obtain 
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[y  - ^(s)]  • [K(s)  - t(s)  Rk  = Rk 

Here  R is  the  radius  of  curvature  of  the  diffracting  edge,  K(s) 
the  unit  normal  vector  at  the  edge,  and  t(s)  the  unit  tangent 
vector  to  the  edge.  Equation  (42)  can  be  written  in  the  form: 

[y  - x(s)]  • t(s)  = [y  - x(s)|cos  6(s). 

The  angle  6 between  the  diffracted  ray  and  h is  given  by 

cos « = . s 

|y  - 5(s) I 

The  radius  of  curvature  of  the  diffracted  ray  R has  been 

a 

defined  above  as  the  distance  |y  - x|  from  the  edge  to  the 
caustic  where  two  adjacent  diffraction  cones  intersect  . This 
distance  is  obtained  by  substituting  Eqs . (44)  and  (45)  into 
Eq.  (43): 

. 2 

R^[cos  6 + P, sin  8]  = Rj^,  sin  8. 

The  result  [Eq.  (46)]  of  the  computation  for  the  radius  of 
curvature  (the  divergence  of  the  elementary  rays),  can  be 
written  in  the  following  form: 

(±)  R,.  sin^B 

R = R j , 

“ cos  6 + Rj^  SsinS 

where  again  Rj^(s)  is  the  radius  of  curvature  of  the  edge,  6(s) 
the  angle  between  the  incident  ray  p and  the  (positive) 
tangent  t(s)  to  the  edge  (i.e.,  the  semi-angle  of  the  diffraction 
: cone)  and  0 • 30/3s  the  angle  rate  between  the  diffracted  ray  and 

I the  edge.  The  unit  vector  Ir  then  is  the  normal  pointing  from  the 

i edge  point  towards  the  center  of  curv'aturc  of  the  edge  in  the 

: plane  of  the  screen.  The  direction  in  which  the  arc  length 

; increases  is  of  no  consequence,  since  both  0 and  sin  0 change 


(44) 


(45) 


(46) 


(47) 
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sign  when  this  direction  is  reversed. 

The  sign  of  the  radius  of  curvature  has  still  to  bo 
determined.  It  is  standard  practice  to  consider  Rj^  as  a basically 
positive  quantity.  If  the  edge  is  curved  away  from  the  source 
point,  one  would  c.xpect  diffraction  to  lead  to  a caustic  focus 
whose  position  is  away  from  the  edge  in  the  direction  of  tlie 
diffracted  rays.  The  wave  should,  therefore,  appear  to  diverge 
from  (r  - ^ ) and  we  have  to  use  the  negative  sign  in  the  above 
formula.  That  this  sign  is  indeed  the  right  one  to  use  is 
borne  out  by  the  results.  The  wrong  sign  would  lead  to  absurd 
results.  We  next  evaluate  the  inverse  1/*^^',  from  Eq . (47) 

, (cos  6 + R.-B  sin  Bj 

1 = - ^ (48) 

a Rj.sin‘‘6 

. " 

We  still  have  to  compute  B.  W'e  have,  from  the  vector  diagram, 
cos  3 = ^ = cos  (p,3s)  (49) 


differentiating  this  with  respect  to  s yields 


- B sin  B 


do 


d"p 
ds  ^ 
2 


(59) 


The  derivative  — y is  computed  with  the  aid  of  the  well 
ds*- 

known  Fresnet  formula  (see  any  text  book  on  differential  geometry) 


If  t is  the  unit  tangent  hector  along  ds  , 

cos  (p  , ds) 


dp 

Ji- 


dp 


= 1, 


cos  B = — 


t, 


and 


dt 


I 


f 


I 
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We  thus  have: 

- 6 s in  6 = 


d^D  d 

57'  ^ 


P • t 


= - ^ (P  • t) 


. 1 


dp 

a? 


1 


-*•  dt 

p • ai- 


p 

1 •?  1 

= p sin  B + ^ cos  (p  ,h)  , 
K 


If  Eqs . (49)  and  (51)  are  substituted  in  Eq . (48),  the  following 
expression  results  for  the  principal  radius  of  curvature 


cos  5 + 1 + 1_  cos  (p  ,ii~) 

Rj^sin^B  ^ sin’B 


1 

’5 

Rj^sin“3 


[- 


cos  6 + 


cos  (p 


In  the  Rubinowicz  notation,  6 = tt  (r,li)  (r  points  from  the 

observation  point  to  edge  element)  . Hence  w'e  have 

i + r-  [cos(r,li)  + cos  (p,ii)] 

a ^ Rj^  sin  B 

The  radius  of  curvature  R^  represents  the  curvature  of  the 
diffracted  wave  front  in  a plane  that  contains  the  edge  element. 
The  last  result  shows  that  diffraction  increases  the  curvature 
of  the  incident  wave  by  an  amount  that  is  proportional  to  the 
curvature  of  the  diffracting  edge.  The  second  principal 
curvature  of  the  wave  front  then  is  simply  given  by  the  inverse 
of  the  distance  from  the  diffracting  edge. 

If  we  enter  this  result  in  Eq . (41i  writing  r (distance 
from  stationary  phase  edge  point)  for  i,  and  recalling  that 
A^/p  is  equal  to  and  sin  B = sin  (p,ds)  we  obtain 


Sir- 


i 

i 

I 


1ms 

I 


(51) 


h]  (5 


(53) 
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, , , . 

1 

cos(r,K)  cos(p,h) 

sin^"  (p,3s)  j 

(54) 


The  Keller  result  is  thus  identical  with  the  Rubinowicz-Kirchoff 
distance  variation  for  a stationary  phase  point.  This  result  is 
quite  suprising.  Even  the  curvature  of  the  incident  wave  front 
has  no  other  effect  than  those  described  by  geometrical  optics 
and  Formats ' principle.  The  Rubinowicz-Kirchoff  does  furnish 
the  directivity  factor 


s in (p ,ds) 


cos (r ,n) 

1 + cos(p,r) 


(55) 


for  the  incident  wave,  and  a similar  factor  for  its  mirror  image, 
i.e.  for  the  reflected  wave.  Since  the  sum  of  the  two  factors 
docs  not  reduce  to  a simple  expression,  it  is  expedient  to  stop 
at  this  point  and  to  consider  the  fields  generated  by  the  incident 
and  the  reflected  waves  separately. 

Thus,  we  have  obtained  the  correct  distance  dependence 
for  the  diffraction  field  in  the  high  frequency  limit  from 
geometric  optical  principles.  But  the  method  gives  no  indication 
of  the  angular  variation  of  the  diffraction  field.  Keller 
determines  the  directivity  factor  by  replacing  the  constant 
with  whatever  is  needed  to  make  his  solution  agree  with  the 
exact  solution  for  the  straight  edge  and  using  the  resulting 
formula  for  all  other  curved  edges.  In  contrast  to  Keller's 
methods,  the  Rubinowicz-Kirchoff  integral  can  be  applied  to 
all  but  the  very  low  frequencies,  i.e.  whenever  the  linear 
dimensions  of  the  diffractor  exceed  one  wave  length.  Likewise 
the  Rubinowicz  integrals  furnishes  the  directivity  function,  or 
at  least  a very  good  approximation  to  it.  It  seems,  however, 
that  the  Rubinowicz-Kirchoff  method  in  present  form  fails  for 
blunt  three-dimensional  diffractors,  which  can  be  handled  very 
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b.  CAUSTICS 

Rays  diffracted  by  adjacent  elements  of  a curved  edge, 
generate  diffraction  cones  which  intersect  along  a second 
order  curve  common  to  both  cones,  as  seen  in  Figure  fll). 

Since  the  wave  length  of  the  incident  field  has  been  assumed 
as  being  infinitely  small  and  since  all  the  contribution  from 
the  edge  element  pass  through  this  curve  which  is  a space  of 
zero  volume,  the  energy  desnity  there  becomes  infinitely  great. 
This  curve  acts  as  if  it  were  a line  focus.  As  we  move  smoothly 
along  the  edge,  these  focal  lines  form  a surface  of  particular 
concentration  of  the  diffracted  energy,  which  we  called  a 
caustic  surface.  The  diffracted  rays  appear  to  originate 
from  this  caustic  surface. 

As  an  example,  we  will  consider  the  focusing  properties 
of  a pencil  of  rays.  All  the  rays  in  a homogeneous  medium 
are  straight  lines  and  are  orthogonal  to  the  wavefronts. 

Consider  now  the  small  pencil  of  rays  which  intersects  a 
wavefront  in  the  small  rectangle  PQRS  (Fig.  12)  bounded  by 
lines  of  curvature.  The  normals  at  F , Q intersect  at  X and 
those  at  R,  S intersect  at  Y.  Similarly  the  normals  at  P,  S 
intersect  at  Z and  those  at  Q , R intersect  at  K.  All  the  rays 
through  the  rectangle,  to  the  first  order  of  approximation, 
pass  through  the  lines  XY  and  ZW  which  are  called  the  "focal 
lines"  of  the  pencil.  This  concentration  of  rays  is  easily 
observed  with  a lens  held  obliquely  in  the  sun.  A small  fan- 
like pencil  cutting  the  surface  along  PS  at  Z.  A small  fan- 
like pencil  in  any  other  direction  is  not  brought  into  focus 
at  all.  The  planes  PXY  and  PZW  are  perpendicular  but  XY  is 
not  in  general  perpendicular  to  PX  nor  ZW  to  PZ . 

Only  in  the  special  case  when  the  wavefront  is  spherical 
do  all  the  rays  of  the  pencil  intersect  in  a common  point;  the 
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pencil  is  then  said  to  be  homocentric  point  source  but  the 
pencil  does  not  remain  homocentric  in  general  after  reflection 
and  refraction.  Instead,  the  radiation  becomes  concentrated 
in  focal  lines. 

If  we  carry  out  the  preceding  construction  for  all  pencils 
of  a wavefront  the  points  Z will  generate  a surface  and  so 
will  the  points  X.  Those  two  surfaces  are  known  as  the  focal 
surfaces  or  caustics.  Every  normal  to  a wavefront  is  tangent 
to  both  the  caustics.  The  caustics  are  the  loci  of  the  two 
principal  centers  of  curvature.  Therefore,  for  a wavefront 
which  is  a surface  of  revolution,  the  caustics  consist  of  a 
portion  of  the  axis  of  rotation  and  a surface  of  revolution 
whose  meridian  section  is  the  evolute  of  the  meridian  section 
of  the  wavefront. 

The  caustic  surface  can  be  defined  by  R = r,  so  that 
1/r  - l/Rgj  = 0 where  is  the  distance  of  the  caustic  along 

the  diffracted  ray  from  the  edge  element.  For  a field  point 
at  a caustic  surface. = 0,  the  standard  stationary  phase 
solution  therefore  becomes  infinite.  The  reason  is  obvious. 

The  Taylor  development  in  the  exponent  is  the  integrad  and 
has  been  broken  off  with  the  square  terms.  Because  4>"  = 0, 
recalling  that  also  t ' = 0 for  a stationary  phase  point,  the 
exponent  is  constant  and  extending  the  limits  to  + «>  and  - “ 
makes  the  integrad  infinite.  For  a simple  caustic  surface, 

4)"  = 0 but  4)"'  ^ 0;  in  fact,  we  may  consider  this  relation 
as  the  definition  of  a simple  caustic.  Near  this  caustic 
4i"  will  be  small,  and  the  Taylor  development  of  the  exponent 
about  the  stationary  phase  integrad  must  be  carried  up  to 
third  order  terms: 


K 


(s  - 


(s  - 


(s 


f 


(56) 


Here,  ' (s^)  - 0,  because  by  definition,  the  stationary  phase 
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point  is  at  s^.  But  because  (s  - is  negative  for  s<s^, 

there  is  another  point  for  which  the  above  development  vanished. 
This  point  is  given  by 


4.’(Sj)  = 0 = (Sj  - i (Sj  - (57) 

or  sj  - Sq  = - ' 

provided  4i"/4>'"  is  small  (say  compared  to  tr)  , the  above  series 
assumes  large  values  for  (s  - s^)>>tt  and  the  limits  of 
integration  can  again  be  replaced  by  ±«>.  The  above  condition 
means  that  the  field  point  should  not  be  very  far  away  from 
the  caustic  surface.  The  stationary  phase  integral  can  then 
be  evaluated  (see  Foundations  of  Acoustics,  page  64).  The 
result  is 


♦ = 


S(s)  e-J''*'")  ds  - 


■ 2 (s  ) 


♦ (S^)  * 4>"(s^)^/3*'”  (s^)2 


00 

cos 

Tffv^-mvj 

2 

0 

. } 

dv 


(58) 


where 


m 


3/3k 


7T  (4> 


;r,y4/i 


(59) 


The  integral  is  the  Airy  integral  (see  loc.  cit)  . To  obtain 
4>"',  we  start  out  with  the  expression 

^ - cos  (;,dl)  = i (p  . t)  (60) 
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Hence  we  obtain 


d'o  . 


ds‘ 


1 

T 

p 


t)  cos  (p,ds)  + i + i 


dt 

ars 


(61) 


and  differentiate  again  with  respect  to  s,  making  use  of  the 
following  relations. 


dp  ^ ^ . 7:  _ :^2 

3?  t,  t • t - t 


1,  t 


dt 

3? 


= 0 


(62) 


^ = cos  (p,ds)  = — *t,  ^ ^ unit  vector  along  main  normal) 

K 


and 


dh  ^ (63) 

3?  ■ T 

where  B is  the  unit  vector  along  the  binormal  and  T the  torsion. 
Also  cos  (t,1s)  = - cos  (p,Ss)  and  sin  (r,^s)  = sin  (p,ds) 


(Eq.  27). 


A similar  expression  is 


obtained 


for 

ds-^ 


Hence 


+ 


. 2 

3 sin  3 cos  3 


<jrp—  [cos  (r  ,5) 
K 


1 

*1 

1 

^ 3 cos  3 

cos  (r,ii)  cf)S  (p,K) 

T 

_r 

r p 

dR 

+ cos  (p,^)]  - [cos  (r,ii)  + cos  (p,K) 


In  deriving  Eq.  (64) , we  have  used  the  stationary  phase  condition 
that  cos  3 = cos  (p,ds)  = - cos  (r,ds). 

If  the  incident  wave  is  plane  and  if  part  of  the  edge  is  an 
arc  of  the  circle  of  radius  a whose  plane  is  normal  to  the 
incident  radiation,  thenp"  and  all  higher  derivatures  of  4>  are 
zero.  The  exponent  then  is  constant  for  the  caustic  point  and 
the  integral  is  equal  to  integral  of  S(s)  over  the  arc  of  the 
circle.  Since  S(s)  is  a slowly  varying  function  of  s , it  can 
also  be  considered  to  be  constant  over  the  range  As  and 


(64) 

] 
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|(()^|  = S(s)  As 

where  As  = a(<{>2  ' ^2  ~ *^1  angle  the  circular  arc 

subtends  with  its  center  of  curvature.  The  quantity  is  the 
distance  on  the  field  point  from  the  arc.  For  points  near  the 
caustic  (near  the  center  of  curvature)  which, by  rotating  the 
di f f ractor,  we  assume  to  lie  on  the  x-axis,  the  projection  of 
the  vector  from  the  field  point  to  the  edge  element  on  to  the 
plane  of  the  diffractor  is  given  by  (see  Fig.  13) 


7 7 7 7 7 7 

r'  = (a  - X cos  (p)  + x^  sin^4i  = a*^  + x^ 


2 ax  cos 


= a^(l  - cos  4>) 

r ' = a (1  - ^ cos  (p)  + ... 

a 

where  x is  the  distance  of  the  projection  of  the  field  point 
from  the  center  of  the  circular  arc  on  the  plane  of  the  circular 
arc,  and  0 the  angular  coordinate  of  the  edge  element  (see 
Fig.  13).  The  distance  of  the  field  point  from  the  edge  element 
is  then  given  by 


+ (a  - X cos  ) 


z*^  + a - 2 ax  cos  + 


where 


2ax  cos  41 
1 “ ^ ~ 


: r (1  - + 

° ^0 


2 2^2 
r * z + a 
0 


because  x<<r^.  The  diffraction  integral  now  simplifies  to 
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S(s)  exp»[jkr^  (1  - ax  cos  4i/r^)]  ds , 


= S(^) 


J(k  ax  cos  <t)/r^ 


(68) 


If  the  curve  is  a complete  circle,  = 0,  = 2it  and 

* = S((^)  e'^^’^o  2iTa 

If  (t>2  - 4>2^  <<  1T , then  cos  4)  = 1 , and 

♦ = S(4»)  • e ^^^o  As 


kxa 

r„ 


(69) 


(70) 


where  As  is  the  length  of  the  arc,  whose  curvature  is  constant 
and  whose  plane  is  normal  to  that  of  the  incident  radiation. 

Keller's  method  of  taking  account  of  the  caustic  can  only 
be  applied  to  the  case  when  the  diffractor  is  a circular  piston. 
He  multiplies  his  stationary  phase  solution  for  the  circular 
aperture  by  a factor  that  is  good  for  the  circular  aperture, 
and,  hence  his  account  for  the  caustic  is  only  good  for  this 
aperture . 


7.  THE  FIELD  AT  AND  NEAR  THE  SHADOW  BOUNDARY  AND  THE  PHYSICAL 

SIGNIFICANCE  OF' EDGE  WAVES 


In  the  Kirchof f -Rubinovicz  theory,  the  diffracted  field  is 
expressed  as  the  cumulation  of  waves  that  originate  at  every 
edge  element  and  is  given  by  Eq.  (21), 
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1 1 cos  (n,r) 

sin  (p,ds) 

1 + cos  (r,p) 


•jk(p  + r) 


ds 


(74) 


edge 

The  factor  sin  (p,ds)ds  measures  the  incident  intensity,  i.e.  width 
of  light  beam  normal  to  direction  of  propagation  that  impinges 
on  ds  . At  the  shadow  boundary,  ^ (n,r)  = Tr/2,Hr,p)  = tr,  so 
that  cos  (n,r)  = o and  1 + cos  (r,p)  = 0 and  the  limiting  value 
for  cos  (n,r)/(l  + cos  (r,p))  turns  out  to  be  infinite.  Physically, 
amplitudes  cannot  become  infinite,  and  therefore  one  suspects 
that  the  Rubinowicz  edge  wave  is  only  a mathematical  fiction 
without  physical  reality.  To  see  what  happens  near  the  shadow 
boundary,  let  us  analyze  the  factor  (see  Fig. 14  ). 


cos  (n,r)  _ r cos  (n,r) 

1 cos  (r,p)  r + r cos  (r,p) 

Let  z be  the  distance  between  the  projection  of  the  field  point 
P'  on  the  edge  and  the  diffracting  edge  element  ds  , y be  the 
distance  between  P and  its  projection  P"  on  the  shadow 
boundary,  and  I the  perpendicular  distance  of  P"  to  the  edge 
(see  Fig.  14).  We  then  have 


2 „2  2 2 

T = SL  + y z 


cos(r.p)  = - ^ 

(76) 

cos(n,r)  = ^ 

Hence  the  factor  becomes, 

r cos  (n,r) _ y 

r - it 

r r cos  (r,P) 

Its  main  contribution  thus  originates  near  r = or  y and  z <<  1. 


May  14,  1973 
EJS:SIH:ADS:lins 


> 


- 43  - 

If  y and  z are  small,  we  have 


/^y"  + 


— T V 

z “ + 


- £ 


, 2 2 
1 y + z 

I Z 


(77) 


We  have  developed  the  square  root  into  a Taylor  series  and 
retained  only  the  low'est  order  terms.  For  z = 0,  and  y = 0 
the  result  diverges.  It  is  important  to  note  that  the  edge 
does  not  consist  of  points  but  of  infinitely  small  elements 
ds  = dz  > 0 that  have  a length  greater  than  zero.  If  we 
integrate  the  edge  wave  for  a very  small  edge  element  Az, 
we  obtain 


^ ^ sin  (p  ,ds)  2£y 


Az/2 


dz 


-r  e 


-jk(r  + p) 


77  ^ (P.ds)  2£ 


y + z‘ 

-Az/2 

T Az/2 


tan'^  - 


•jk(r  + p) 


yJ-Az/2 


y sin  a.dt)  e-f'C-  • - 


e-jk(r  * P)sin(p,<is)  (78) 


because  y -►  0 at  the  shadow  boundary,  and  Az  is  very  small  but 
finite  and  also  z has  been  assumed  to  be  small.  Also,  because 
the  integration  interval  Az  is  very  small,  the  exponent  can  be 
assumed  to  be  constant  during  the  integration. 

The  main  contribution  to  the  field  at  the  shadow  boundary 
is  thus  generated  by  the  edge  element  that  casts  its  shadow 
to  the  field  point  at  the  shadow  boundary.  The  integrand  be- 
comes infinite  for  their  midpoint  of  this  element,  but  the 
average  of  i (like  the  integral  over  a delta  function)  is 
finite.  The  contribution  of  this  element  then  is  equal  to 


* 
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half  the  amplitude  of  the  incident  wave  at  the  edge  point.  If 
we  had  applied  the  stationary  phase  method  in  a form  that  is 
valid  when  the  integrad  has  a pole  near  the  stationary  field 
point,  wo  would  have  integrated  over  an  edge  length  of  many  wave- 
lengths, so  that  phase  cancellations  w’ould  occur,  and  would 
have  obtained  the  same  result  for  the  stationary  phase  con- 
tribution except  that  the  distance  factor  ^ would  hav^e  been 

replaced  by  ^ ^ ^ , so  that  the  diffracted  amplitude  is  exactly 
half  the  amplitude  that  is  incident  at  the  field  point  and  not 
half  of  the  amplitude  that  is  incident  at  the  edge,  (see  E.  J. 
Skudrzyk,  1971,  page  528).  Thus, approaching  the  shadow  boundary 
from  the  insonified  zone  leads  to  an  edge -generated  wave  whose 
value  is  exactly  equal  to  half  the  incident  amplitude  and  op- 
posite in  sign  to  it. 

The  factor  of  the  integrand 


1 

r [1  + cos  (r,p)] 


1 

r + r cos  (r,p) 


(79) 


is  the  distance  from  the  edge  to  the  field  point  minus  the 
projection  of  this  distance  on  the  shadow  boundary.  One  is 
tempted  to  interpret  r'  as  the  effective  distance  the  edge 
generated  wave  has  travelled.  For  a point  at  the  shadov\r 
boundary,  ( ^ (r,o)  = tt),  r'  = 0,  for  points  away  from  the 
shadow  boundary  in  the  space  opposite  the  sound  source  r'  < r, 
and  for  points  in  the  insonified  side  of  the  screen  r'  > r. 
However  r'  is  not  the  distance  from  any  particular  point  on 
the  shadow  boundary.  Also,  the  location  for  the  "source  point" 
at  r'  » 0 turns  out  to  be  a function  of  the  field  point,  r, 
even  if  we  move  along  a line  passing  through  this  "source  point" 
and  the  diffracting  edge. 

Thus  it  is  not  possible,  therefore,  to  interpret  the 
diffracted  field  as  resulting  from  elementary  waves  that 
originate  at  the  edge.  The  phase  of  the  diffracted  field 
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does  agree  with  that  for  an  edge -generated  wave,  but  the 
amplitude  of  the  diffracted  field  varies  in  a strange  manner 
with  the  distance  from  the  edge.  The  elementary  edge  waves 
in  the  Rubinowicz  integrand  thus  represent  only  a mathematical 
formalism  with  no  physical  significance.  The  elements  of  a 
diffracting  edge  do  not  generate  real  waves  that  propagate 
from  the  edge  points  as  spherical  waves  to  the  field  point 
like  the  waves  in  Huyghens  ' principle.  However,  the  Rubinowicz 
integral  proves  that  each  element  contributes  to  the  diffracted 
field  observed  at  the  field  point,  and  that  the  phase  of  this 
field  is  the  same  as  that  of  a wave  that  travels  from  the 
source  to  the  diffracting  edge  element,  and  from  there  to  the 
field  point.  However,  if  we  add  up  the  contributions  of  all 
the  edge  elements,  real  waves  result  that  originate  at  the 
stationary  phase. 

According  to  the  transformed  Kirchoff  integral,  the  shadow 
boundary  is  the  exact  set  of  the  sources  that  give  rise  to 
the  scattered  field.  The  contribution  of  the  shadow  surface 
can  then  be  transformed  into  that  of  a source  distribution 
around  the  edges  of  the  diffractor.  However,  these  sources 
are  of  a non-physical  nature  and  have  properties  that  cannot 
be  realized  by  real  sources.  For  example,  there  are  no  known 
sources  nor  any  simple  multiple  sources  whose  contribution 
varies  with  the  distance  as  the  integrand  varies  in  the 
Rubinowicz  integral. 

This  result  explains  why  the  theory  of  diffraction  based 
on  equivalent  elementary  edge  sources  has  not  been  more 
successful  in  general,  and  why  it  had  not  been  possible  to 
represent  most  of  the  exact  solution  in  a similar  form.  The 
edge  sources  are  only  a formal izm  equivalent  to  the  Kirchoff 
theory,  and  equivalent  to  the  exact  solution  for  straight  edges 
and  wedges . 
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8 . THE  STR.MGHT  EDGE  - KIRCHOFF  AND  SONIMERFELn 

For  a plane  wave  incident  at  an  angle  with  the  plane 
that  contains  the  edge  whose  wave  fronts  are  parallel  to  the 
edge,  the  stationary  phase  point,  (see  Fig.  15),  the  field 
point  and  incident  ray  are  in  the  same  plane,  ® and 

^ (r  ,p)  = 360°  - ((J)  - <t>^)  , 

^ (n  ,r)  = <{i  - 4)p  - 90°  , 

sin  (o,ds)  = 1 


The  angle  factor  thus  reduces  to 


cos  (n,r)  sin  ro,ds) 
1 + cos  (r,p) 


sin  (4)  - 4)^) 


and  the  complete  solution  is 

' -j  (kr  + tt/4) 
$ = - £ 

2/^Trkr 


Here  {i4  - 4>q)  represents  the  angle  of  the  rays  from  the  edge  to 
the  field  point  relative  to  the  rays  that  delineate  the  shadow 
boundary  of  the  incident  wave.  The  exact  Sommerfeld  solution 
leads  to  the  same  first  factor,  but  the  second  factor  is 

1/cos  [ (4»  - 4>^)/2].  similar  result  then  is  obtained  for  the 
diffraction  generated  by  the  reflected  wave.  It  is  represented 
by  a similar  solution  except  mot  the  angular  distance  from  the 
shadow  boundary  - 4>  now  is  replaced  by  that  for  the  reflected 
wave  4)  4>q  - 2tt  (sec  Fig.  15).  The  angle  (4  + 4>q  - 2tt)  then 

represents  the  angle  of  the  field  point  relative  to  the  shadow 
boundary  (the  angular  limit)  of  the  reflected  wave. 

The  Kirchoff  formula  gives  the  same  distance  dcpcndancc  as 


(80) 

(81) 

(82) 

(83) 

(84) 
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the  exact  theory,  but  the  directivity  function  is  in  error  by 
the  factor  sin  [(:>  - This  error  is  less  than  30  percent 

if  the  angle  betAvecn  field  point  and  shadow  boundary  is  less 
than  90°.  Since  the  diffraction  at  greater  angles  is  dominated 
by  that  of  the  reflected  waves  which  is  well  approximated  by 
the  Kirchoff  solution  in  this  angular  range,  the  Kirchoff  theory 
may  be  considered  as  a very  good  approximation  to  the  exact 
theory  of  diffraction.  However,  it  is  necessary  to  include 
the  reflected  wave  field  and  its  diffraction  separately. 

If  the  incident  wave  front  is  not  parallel  to  tlie  edge, 
so  that  the  incident  rays  make  an  angle  6 = ^ (p,ds)  with  the 
edge,  then  the  stationary  phase  point  is  no  longer  the  point 
on  the  edge  that  is  closest  to  the  field  point,  but  is  displaced 
from  this  point  towards  the  direction  of  the  source.  To 

find  this  point,  we  begin  by  drawing  the  line  from  field  point 

that  perpendicularly  intersects  with  the  edge.  This  point  P 
will  be  designated  as  the  origin  of  s (see  Tig.  16). 

The  incident  ray  that  corresponds  to  this  origin,  s = 0, 
is  called  o^.  The  distance  from  the  source  to  the  edge  plus 

that  from  the  edge  to  <^ield  point  is  a function  of  s and  can 

be  written  as  ' 


p(s)  + r(s)  = Pq  - s cos  6 


To  find  the  specific  location  of  point  (sj^  p for  which  the 
phase  is  stationary,  we  set 


'•j  = 0 - - cos  8 • 

dS 


/FT 


"T 


= 0 


Hence 


cos^B  (d^  + s^)  = s^  , 


(85) 


(86) 
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d cos  B = (+)  s sin  6 , 


or 


tan  6 = d/(s)g  p 


d cos  3 
s in  B ’ 


(87) 


as  would  have  followed  directly  from  the  edge  law.  The  quantity 
(s)  represents  the  distance  of  the  stationary  phase  point 

from  the  point  on  the  edge  that  is  closest  to  the  field  point. 
Note  that  the  stationary  phase  point  is  very  far  away  from  the 
point  at  the  edge  that  is  closest  to  the  field  point  unless  the 
incident  rays  are  normal  to  the  edge.  For  instance,  if  the 
incident  rays  form  an  angle  of  45°  with  the  edge,  the  stationary 
phase  point  is  as  far  away  from  the  point  on  the  edge  that  is 
closest  to  the  field  point  as  the  field  point  is  from  the  edge. 
Hence,  for  the  stationary  phase  point 


(P  ^ r)^ 


.p. 


d cos^B  ^ 
sin  B 


d 


cos  6 
T~ 

sin“B 


1/2 


= P, 


s in  B 


(1  - cos^B) 


+ d sin  B 


(88) 


The  right  hand  side  represents  the  exact  distance  from  the 
source  (very  far  away)  to  the  stationary  phase  point  and  from 
there  to  the  field  point.  Thus  we  obtain  the  result: 

$ = - B g-ik  (p^  + d/sin  B)  + Jtt/4 


<t>  - <P, 


(89) 


where  r = d/sin  6,  and  d is  the  perpendicular  distance  of  the 
field  point  from  the  edge.  represents  the  phase  of  the  incident 

wa%re  at  the  origin  s = 0,  which  is  the  point  on  the  edge  defined 
by  the  projection  of  the  field  point  on  the  edge. 

If  the  incident  wave  is  a spherical  wave,  if  s^  is  the 
s coordinate  of  the  normal  projection  of  the  source  point  on 
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the  edge,  and  is  the  distance  from  the  source  to  the  projection 
point  or  the  edge,  then  the  s cordinatc  of  the  stationary 
phase  point  satisfies  the  relation  (see  fig.  17) 


and 


tan  6 = ^ = 


s s - s 
0 


(S) 


s . p , 


^ ^0 

p + J 
o 


(p  + r) 


"o  - = 


s .p, 


cos  ^ cos  0' 


cos 


(Po  - ^ s^ 


-n  1/2 


= R 


(90) 


where  R is  an  abbreviation  for  the  sum  [p  + r] 


s.p 


We  now  have 

-j  (kR+r/4)  ^ 

♦ j sin  6 tan  — (91a) 

2/2 irk  pr 


where  p,r  represent  the  distance  vectors  to  the  stationary 
phase  point,  and  p,r  the  corresponding  distances.  Thus,  the 
total  Kirch off  solution  for  the  diffracted  field  by  the  incident 
and  reflected  fields  give: 

- j (kR  + tt/4) 

<!>,=  - sin  S 

^ 2/TnY'pTT 


4>  - 


tan 


tan 


(91b) 


The  directivity  factors  of  the  diffracted  field  in  the  Kirchoff 
solution,  equation  (91b)  for  rigid  (+ j and  pressure  release  (-) 
edges,  are  shown  in  Figures  (18)  and  (19)  respectively  for  a few 
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incidence  angles.  The  corresponding  total,  exact  (Sommerfeld) 
solution  for  the  diffracted  field,  see  Bowman,  Senior  and 
Uslenghi,  1969,  can  be  written  as: 


4 


d 


- j (kR  + tt/4) 

sin 

2/2Trk  pr  k" 


1 


cos 


4>-<l 

~T 


1 

cos  — 


(91c) 


The  directivity  factors  of  the  diffracted  field  in  the  Sommerfeld 
solution,  equation  (91c)  for  rigid  (+)  and  pressure  source  release 
(-)  edges,  are  shown  in  Figures  (20)  and  (21)  respectively. 

The  backscattered  field  when  3 = tt/2,  r = p = d = , 

Sq  = o,  R = 2p,  and  4>  = 4>q  becomes 


f + tan  (Kirchoff) 


'bs 


• j (2kp  + tt/ 4) 

4p/rn<p 


(91d) 


1± 


cos  4> 


(Sommerfeld) 


The  directivity  factors  of  the  backscattered  field  for  the  Kirchoff 
and  Sommerfeld  solutions  are  shown  in  Figure  (22)  for  rigid  edges 
and  in  Figure  (23)  for  pressure  release  edges. 


9.  THE  WEDGE 


Since  the  Kirchoff  solution  depends  only  on  the  shadow 
boundary,  it  is  the  same  for  both  a wedge  and  a straight  edge. 
For  a plane  incident  wave,  the  diffracted  field  of  a wedge 
becomes : 
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^Kirchof f 


Z/TtHTt 


<fc  - (j) 


cos 


(92) 


In  the  exact  solution,  tlie  factor  in  the  braces  is  replaced  by  unity. 

A wedge  is  characterized  by  the  apex  angle  2C1  of  the  wedge,  and 
the  angle  vn:  that  supplements  it  to  Zv:  is: 


V7T  = 2tt  - 2Q  , 

V7T  being  the  number  of  radiants  exterior  to  the  wedge. 

For  a straight  edge,  = 0 and  v = 2. 

For  a wedge  of  angle  90°  , Q = tt/4  and  vtt  = 2tt  - tt/2  = ^ , v = j , 
likewise  for  wedge  of  180°,  SI  = ■n'/ 2,  vtt  = 2tt  - tt  = tt  , and  v = 1 (93) 


For  a wedge,  the  exact  solution  for  the  diffraction  field  of  an 
incident  plane  wave  is  given  as  (J . J.  Bowman,  I.B..A..  Senior  and 
P.L.F.  Uslenghi  , 1969)  : 


4>o 


j (l<r  + tt/4) 


2 /2  TT  k r 


2 . TT 

— Sin  - . 
V vj 


cos  - - 


cos 


1 

TT 

cos  — + cos 

I 

V 

L ^ y J 

= t wCT„f)  (94) 

where  = '^*'^0  angle  the  observer  makes  with  the  shadow 

boundary  of  the  incident  wave  and  angle  the 

observer  makes  with  the  shadow  boundary  of  the  imago  wave,  see 
Figure  (24)  . 

The  directivity  function  of  the  diffracted  field  from  a wedge, 
given  in  the  bracketed  terms  of  equation  (94)  , is  plotted  for 
wedge  angles  ZQ  = 30°,  60°,  90°,  and  120°  for  rigid  (+)  in  Figures  (25, 
27,  29  and  31),  and  for  pressure  release  (-)  in  Figures  (26,  28,  30, 
and  32),  respectively. 

The  backscattered  field  of  these  wedges,  when  <p  = 4i^,  were  plotted 
in  Figures  (33)  and  (34)  for  rigid  (+)  and  pressure  release  (-) 
respectively.  It  should  be  noted  that  the  diffracted  field  represented 
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by  equation  (94)  becomes  unbounded  at  the  following  angles: 


'V  = = 71 

4,-4, 


(j,  + ^o  = 3--2n  I foj.  2J1  ,<  7T 


(j)+  4 = 7T  + 

o 


provided  that 
n 4 41  ■<  360° 


The  factor  ^ sin  — , shown  in  Figure  (35),  can  be  interpreted 
as  the  discontinuity  factor  for  a wedge.  For  a blunt  wedge 
(straight  edge),  the  factor  is  unity.  The  factor  for  an  almost 
flat  wedge,  whose  angle  is  20  = tt  - e , or  whose  angle  v is: 


V = 2.?_.:__2g  = 1 ^ L 

7T  IT 


become  for  small  e , 


— sin  — = 2(1  - — + ...)  sin  (tt  e)  = 2e  ... 

V V ^ TT  ■'  ^ 


(95) 


it  is  proportional  to  the  change  of  the  slope  multiplied  by  two. 


1 0 . COMPARISON  OF  THE  KIRCHOFF  STR.AIGI1T  FDGF.  SOLUTION’  KITH 
THE  SOMMERFELD  STRAIGHT  EDGE  .AND  WEDGE  SOLUTION 


The  solution  for  the  straight  edge  consists  of  the  linear 
combination  of  two  identical  terms  with  different  arguments: 

!refu) 


*d  " 


- 


± w 


w 


± w 


(96) 


where  represents  the  angular  distance  from  the  shadow 

boundary  of  the  incident  wave  and  T n * 2tt  that  from  the 
shadow  boundary  of  the  reflected  (the  image)  wave  (see  Fig.  18  § 24). 
The  plus  sign  applies  for  a rigid  straight  edge,  the  minus  sign 
for  a pressure  release  straight  edge.  The  two  terms  a^'e 
identical  except  for  a difference  in  the  origin  of  the  angle 
count.  It  is  therefore  sufficient  to  consider  only  one  term 
alone.  The  distance  factors  f(r)  are  the  same  in  the  Kirclioff 
solution  and  in  the  Sommerfeld  solution.  The  only  difference 
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is  in  the  angle  functions.  The  angle  function  for  the  Kirchoff 
solution  is  tan  T/Z,  and  that  for  the  Sommerfeld  solution 
is  l/(cos  T/Z}.  Fig.  36  shows  a comparison  of  these  two 
functions.  In  the  regions  where  scattering  is  significant, 
i.e.  near  the  shadow  boundary,  the  two  solutions  are  practically 
identical.  The  two  functions  differ  for  scattering  in  directions 
normal  to  the  shadow  boundary  by  the  factor  sin  45°  = 0.707  or 
30|.  Thus  for  diffraction  into  the  half  space  that  extends 
in  the  direction  of  propagation,  the  Kirchoff  error  is  practically 
negligible.  The  Sommerfeld  value  then  is  1.4,  the  Kirchoff 
value  is  1.  For  the  backscattered  amplitude,  4'=0  and  the 
Kirchoff  backscattered  amplitude  is  zero,  while  the  Sommerfeld 
backscattered  amplitude  is  unity.  Backscattering  is  predominent ly 
determined  by  the  diffraction  field  of  the  reflected  wave.  We 
can  thus  conclude  that  the  Kirchoff  error  for  plane  diffractors 
is  negligible  under  all  practical  situations  provided  the  radius 
of  curvature  is  large  compared  to  the  wave  length. 

The  Kirchoff-Rubincwicz  solution  for  the  straight  edge 
leads  to  the  angle  factor  tan  [(4)  - 4)^)/2],  where  the  exact 
Sommerfeld  solutions  differ  from  that  of  the  Kirchoff  solution 
practically  only  in  the  range  (4  - > fr/4,  or  (J)  - > tt/2, 

i.e.  in  the  range  in  which  scattering  is  extremely  small.  In 
this  range  the  diffracted  field  is  practically  exclusively  that 
of  the  reflected  wave.  We  may  attempt  to  improve  the  solution 
using  Keller's  method,  and  simply  drop  the  factor  sin  (4  - 4>q)/2 
in  the  Kirchoff  solution  in  all  cases  of  edge  diffraction. 

In  the  Kirchof f - Rub inowicz  theory,  the  diffraction  phenomenon 
is  represented  by  an  integral  over  the  shadow  boundary.  The 
resulting  integral  then  is  transformed  into  a line  integral 
over  the  edge  that  generates  the  shadow  boundary.  The  solution 
reproduces  the  correct  distance  dependence  of  the  diffraction 
field  at  high  frequencies,  and  is  dimensionally  correct.  As 
has  been  proven  in  the  literature  that  for  plane  diffractors. 
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the  Kirchoff  solution  represents  a very  good  approximation  in 
the  angular  range  near  the  direction  of  incidence,  provided  the 
linear  dimensions  of  the  diffractor  exceed  one  to  two  wave 
lengths,  (see  Figs.  3 and  4),  We  can  therefore  assume 
that  the  Kirchof f-Rubinowicz  solution  also  represents  a very 
good  approximation  to  the  distance  dependence  of  the  exact 
solution  of  Sommerfeld. 

The  Sommerfeld  computation  does  not  include  the  terms 
that  contain  the  curvature  R,,  of  the  edge.  It  would  be  very 

Iv 

difficult  to  derive  an  exact  solution  for  a curved  edge  at 
low  frequencies  when  the  edge  is  not  very  long  compared  to 
the  wave  length.  We  may  expect  that  the  Kirchoff  solution 
without  the  factor  sin  [ (4>  - 4i^)/2]  will  give  a very  good 
approximation  to  the  exact  solution  for  edges  with  slight 
curvature . 

It  will  be  recalled  that  the  Kirchoff  theory  yield  a similar 
formula  except  for  an  additional  directivity  factor  sin[(ij)  - 4>q)/2] 
for  the  diffracted  field.  A similar  angular  factor  relative  to 
the  shadow  boundary  is  associated  with  the  reflected  wave,  but 
the  Kirchoff  theory  does  not  contain  the  factor  ^ sin  ^ of  the 
exact  Sommerfeld  result.  In  the  exact  theory,  the  angle 
space  outside  the  edge  Zir  - 2^  is  dilated  to  Ztt  by  dividing 


((i  by  V where  2t  - 2J1  = vtt  , 


T . 
me 


in  terms  cos 


'V  . 
me 


- 251 


= 2tt 


so  that  the  argument 


fZTT  - 


and  that  of 


ref  1 


in  the  term 


2tt  - t 


cos 


ref  1 


alxtfays  varies  from  0 to  2v 


(9-/) 


For  wedge  angle  of  90°,  we  see  that: 


251  = 


TT 

I 


f 


2 

— sm 

V 


TT 

V 


4 

I 


s in 


2tt 

3"  ■ 


1.17. 


V 
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“ 

Similarity 

, for  wedge 

angle 

135°  , 

we  have: 

5 2 

V = T » — 

4 V 

TT 

sin  — 

V 

II 

oo 

4 

sin  ^ = 

8 

S """ 

The  deviation  between  Kirchoff  and  exact  theory  become  serious 
when  the  wedge  angle  2f2  exceeds 

11.  DIRECTIVITY  FACTOR  OH  THE  EDGE  RADIATION 


The  diffracted  field  establishes  the  mathematical  continuity 
from  the  illuminated  to  the  shadow  zone  (see  Sec.  6).  The  peak 
of  the  diffracted  field  propagates  in  the  same  direction  as 
the  incident  wave  and  is  half  as  great  as  this  wave. 

Thus  the  diffracted  field  near  the  shadow  zone  becomes 
(see  Eq.  32) 


»d  ■ 


1 * ‘’’sin  e 


ro/  i . i , 

r p RJ 


t(*) 


(98) 


where 


cos  (o,h)  + cos  (r,h) 

y z 

sin  (p,ds) 


and  for  the  Sommerfeld  solution,  the  angular  factor: 


f(«^)  = 


1 


-1 


1 


cos 


e ^ " sin  (e/2) 
o cos  ■ ■ '2‘  * 


(99) 


2 

where  in  the  vicinity  of  the  shadow  boundary 

^ - (^0  = tt  ♦ e (100) 

Here  e is  the  angle  or  latteral  space  deviation  from  the  shadow 
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boundary.  The  above  expression  for  f(4')  applies  everywhere 
with  the  exception  of  the  immediate  vicinity  of  the  shadow 
boundary.  At  the  shadow  boundary,  its  value  is  given  by  the 
amplitude  of  the  diffracted  field  which  is  equal  to  half  that 
of  the  incident  field  (see  Sec.  6).  Hence 


" P'>  sin  6 


rp 


/TTTTx 


f(4>) 


eikR 

nr 


where  R is  the  distance  between  source  and  receiver  and  is 
equal  to  r + p and  p and  r the  distances  of  the  stationary 
phase  point  from  source  and  receiver.  Thus, 


f«>)  = 


R 


1 < 2_ 
sin  e/2  e 


where  e = e is  the  smallest  value  of  e,  we  are  allowed  to 
0 

use  in  Eq.  (1023  - Eor  this  value  of  c,  the  diffracted 
amplitude  is  equal  to  half  tlie  amplitude  of  the  incident 
wave.  The  amplitude  then  stays  at  this  value  and  as  e 
passes  through  the  shadow  boundary  to  - then  decreases 
again  proportionally  to  1/e.  We  thus  have  e > e^ . 

To  simplify  the  computation,  let  the  incident  wave  be 
very  nearly  a plane  wave  (p  -*■  <»)  and  let  it  be  normally 
incident  on  the  edge.  Let  R^,  = ®.  The  stationary  phase 
point  then  is  concident  with  the  projection  of  the  field 
point  onto  ihe  edge.  To  fix  the  amplitude  of  the  incident 
wave,  let  p r.  We  then  have 


(101) 


(102) 
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f(<t)  = 


JlX^  2TTk 


2r 


Akr 


For  values  of  e greater  than  = l/sin(e/2)  = 2/e 

up  to  crudely  e/2  = it  , e = 2tt  , i.e.  for  the  complete  angular 
range.  For  p = r,  we  thus  have 


e > e. 


P 

o 

TF 


p 

o 

p Akr  e 


P 


incident 


1 

Akr  e 


Thus  the  scattered  pressure  decreases  inversely  proportional 
to  the  perpendicular  distance  e from  the  shadow  boundary. 


12.  EDGES  OF  FINITE  LENGTH 


The  theory  of  geometrical  diffraction  is  based  on  Fermat's 
principle,  see  Appendix  III,  which  is  equivalent  to  the  presence 
of  a stationary  phase  point  well  within  the  length  of  the 
diffracting  edge.  The  diffraction  field  is  then  considered 
to  arise  primarily  from  the  contributions  of  the  edge  elements 
at  and  near  the  stationary  phase  point.  If  the  edge  is  not 
long  enough,  or  is  so  situated  relative  to  the  field  point 
that  no  stationary  phase  point  is  located  within  the  length 
of  the  edge,  then  no  diffracted  field  will  reach  the  field 
point . 

For  the  case  when  the  stationary  phase  point  is  within 
the  length  of  the  edge,  the  exponent  function  in  the 

diffraction  integral 

* 

I.  S(s.r)  " »>ds, 

edge 


(103) 


(104) 


(105) 


J 
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can  be  developed  into  a Taylor  series. 

The  stationary  phase  point  satisfies  the  condition 

g 

jg-  (p  + r)  = 0,  which  will  be  considered  as  the  origin 

for  the  s-coordinate  along  the  edge.  The  series  expansion 
of  the  exponent  about  the  stationary  phase  point  s = 0,  where 
r = and  p = becomes : 


(p  + r)  = 


(Po  ’"o^  3^Cp  + r) 


3s' 


P 

O ’ O 


If  the  sound  impinges  normally  to  a straight  edge 


2 2^2 
P = Pq  ^ s , 


P = Pq  ^ J ^ 


2 2,2 
r = r + s 

O ’ 


r - i s^/ro  ♦ 


3 (p  r) 

Vs 


Po’^o 


L_  . L_  . 

Po 


(106) 


(107) 


(108) 


and  the  exponent  in  the  integrand  in  Equation  (105)  above 
reduces  to  the  following  Taylor's  series: 


-jk 


(109) 


Since  the  amplitude  factor 
as  slowly  varying  function 


S(s,r)  of  the  integrand  is  assumed 
of  s,  the  integral  is  basically 
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•7 

Fresnel  integral  of  the  variable  v*" 


(see 


E.  Skudrzyk,  1972,  page  122).  The  value  of  the  integral  in 
Equation  (105)  is  given  by  twice  the  distance  from  the  origin 


of  the  Cornu  spiral  to  the  point  v^  = ^ s 

A. 


max 


tnat 


represents  the  upper  limit  of  integration  of  Equation  (103). 
Except  for  insignificant  fluctuations  due  to  s being  located 
on  the  spiral  near  its  upper  end,  see  Figure  38,  the  integral 


assumes  its  final  value  1 


2 e if  V > /t/2. 


If  the  projector  and  receiver  are  in  the  central  plane 
of  symmetry,  the  solution  for  a finite  edge  of  length  ^ ~ » 

therefore  is  the  same  as  for  the  infinite  edge  only  if  we 
satisfy  by  the  following  limits. 


1 

I 


k 


h-  > TT 

T~  ■ T 


(110) 


so  that  the  finite  limits  of  the  integral  (105)  can  be  replaced 
by  infinite  limits  of  the  length  cf  the  edge. 

For  backscattering , = r^  , the  maximum  limit  range  r^^  for 

the  geometrical  theory  of  diffraction  to  be  valid  is  reached  when 


or 


where  (h/X)  and  (r'/X) 

m 

edge  and  range  normalized  to 


are  the  dimensionless  length  of 
the  acoustic  wavelength.  lor 


(IIF/ 


60 
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example,  in  water  at  30  kHz,  X = ^ m,  and  we  obtain  the 
following  table: 


^m 

r ’ 
m 

T 

h 

Y 

h 

40 

800 

20 

1.0 

1.0x10^ 

Zxio"* 

100 

5.0 

4.0x10^ 

Sxio'^ 

200 

10.0 

1 . 6x10"^ 

3.2x10^ 

400 

20.0 

As  the  size  of  the  target  is  reduced,  the  values  of  h become 
smaller  and  the  exponential  is  practically  equal  to  unity.  The 
integral  in  (105)  then  is  no  longer  a stationary  phase  integral 
but  becomes  proportional  to  s. 

If  the  angular  range  that  contributes  to  the  intensity 
at  the  field  point  P is  26,  (see  Fig.  39)  then  we  hav'e 

* a h X 

® rr:  ■ if 

m 

6 = tan  0 = sin  8 = ^ (11 

The  result  is  similar  to  the  diffraction  formula  for  a slot  or 
an  aperture,  except  for  a slightly  different  numerical  constraint. 

We  may  conclude  that  either  all  diffracted  beams  open  up  to 
an  angle  6,  or  else  there  w'ill  not  be  any  elementary  beam 
contributing  at  the  field  point  if  the  distance  r is  greater 
than  that  given  by  the  above  formula,  Equation  (110). 
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It  is  interesting  to  note  that  physically,  a diffraction 
field  from  finite  edge  acts  similar  to  the  radiation  field 
from  a line  array  of  finite  length;  both  generate  a beam 
pattern  whose  angular  opening  is  described  by  the  above 
formula . 

For  a plane  incident  wave,  let  Pq  ^ in  Equation  (110), 
and  the  limit  range  r^^  of  the  geometrical  theory  of  diffraction 
is  given  by 


Similarly,  the  angle  range  that  contributes  to  the  intensity 
at  the  field  point  P is  6 , see  Figure  (3S)  , now  given  by 


tan  e = ^ 

o 

• • • X 

9 = tan  6 = sin  ^ 


13.  THE  LIMITS  OF  GEOMETRICAL  REFLECTION  AN’D  DIFFRACTION  FOR 

F 1 NTrn : f l at"  r E fI,~f: ctors 

It  is  well  known  that  a plane  rigid  reflector  is  equivalent 
to  a planar  piston  projector.  Every  point  on  the  rigid  surface 
of  a reflector  can  be  considered  to  have  a particle  V'elocity 
that  is  180°  out  of  phase  with  the  incident  wave.  Thus,  if 
Pj  (x,y)  is  the  amplitude  of  an  incident  pressure  wa\e  at  the 
surface  of  a plane  reflector,  then  the  reflected  far-field 
pressure  amplitude  P(r,9,4))  at  a field  point  (r,9,4))  is  given 
by  the  Huygens -Ray leigh  integral: 


(113) 


(114) 
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! 


Cx.y)  c - 

dx  dy  (115) 

r 

where  (x,y)  are  the  coordinates  of  a point  in  the  plane  reflector, 

(dx  dy)  is  an  elemental  area  of  the  plane  reflector,  and  r is 

the  distance  of  the  field  point  to  the  point  (x,y)  on  the  surface,  a, 

of  the  reflector,  sec  Figure  (40)- 

If  the  origin  of  the  coordinate  system  (x,y)  is  located 
on  the  reflector,  whose  distance  from  the  field  point  r^  is 
given  by 


P ( r , 6 , 4i)  - 


'7  2'’'^ 

r“  = x +y‘'  + z", 
o o o o’ 


then  the  distance  r from  the  field  point  to  any  point  (x,y) 
on  the  reflector  becomes  , see  Figure  (40)  , 


">  7 7 7 

r*"  = (x-x  ) + (y-y  ) “ + z 


2 2 2 
= p + r — 2xx  — 2yy  = r 
o 0 ' o o 


o2  2xx^ 

P j.  "I  P 

T ^ i ’I — 


2\'v 


' • 0 


where  p is  the  distance  of  the  point  (x,y)  to  the  origin 


The  direction  cosines  of  the  line  connecting  the  center  of  the 
radiator  with  the  field  point  are 


Xq/To  * cos  (r^,x)  » cos  9,  Y^/r^  = cos  (r^,y)  = cos  4> 


and 


r * r 


o/ 


n -J - 

1 " — 7^  cos  6 

o 


"Tv" 


cos  4> 


(116) 
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If  the  linear  dimcnsiuiis  of  the  reflector  are  small  compared 
to  the  distance  from  the  field  point,  p /r^“  «1,  the  square 
root  can  he  developed  into  ' power  series 


r 


1 - T 


[(x/r^)  cos  0 + (x/Tq)  cos  t]  ♦ 


1 

X cos  0 + y cos  i) 

2 

•f 

* 

I 

r 

o 

J 

. 

Thus,  retaining  terms  up  to  the  second  power  in  the  exponent  of 
the  integral  in  equation  (115),  one  obtains: 


PCTq  ,9  ,<f) 


Ttt  r 


Pj  (x.y) 


exp  jk 


X cos  6 + y cos 


[x  cos  0 


cos 


dx  dy , 


where  r =r  over  the  surface  of  the  reflector, 
o 

Let  an  incident  plane  pressure  wave,  whose 
amplitude  at  the  surface  of  the  reflector  is  a constant  , 
insonify  a rectangular  plane  reflector.  If  the  coordinates  of 
the  reflector  extends  over  -a  < x < a and  -b  < y < b , then  the 
reflected  pressure  amplitude  at  a field  point,  to  a first  order 
(Fraunhofer)  approximation,  is  given  by  the  integral 


P(r  ,0 ,4>) 


jk  P„ 

Tn 


- j kr 
e ■'  o 


« + a 


- a 


fb 


-b 


^jk[x  cos  0 + y cos  ;>] 


jk  (4ab)e  ^^’"o 


D(0,<^)  , 


(117) 


(118) 


dy 


(119) 
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where  directivity  factor  D(9,4))  is  given  by 


D(e,<t)  = sin  (ka  cos_  9) 
ka  cos  0 


sin  C k h cos 
kb  cos  (J 


(120) 


V 


i 


It  is  obvious  that  maximum  reflection  is  attained  at  a field 
point  in  a direction  normal  to  the  plane  reflector  at  the  origin, 
i.e.  as  0 and  ^ , where  tlie  directivity  factor  D(0,4i)  ap- 
proaches unity.  The  reflected  pressure  at  a far  field  point, 
where  r^  >>  a and  r^  >>  b,  decreases  as  l/r^,  just  like  spherical 
spreading  of  a point  source.  Thus,  for  such  large  distances, 
the  laws  of  geometrical  optics,  whose  pressure  amplitude  decreases 
as  l//r^  , are  no  longer  applicable. 

The  proceeding  computions,  see  section  12,  do  not  apply 
for  ranges  r < r^^ , where  r^  is  the  limit  range  of  the  geometric 
theory  as  defined  in  liquation  (113).  Thus,  one  needs  to  get  a 
better  approximation  by  retaining  more  terms  in  the  eikonal 
expression  of  Equation  (118).  For  a field  point,  whose  projection 
on  the  plane  reflector  falls  within  its  boundaries,  let  the 
origin  of  (x,y)  coordinates  be  that  projected  point.  Thus, 

0 = 0 , <f>  = 0 and  the  integral  for  the  reflected  pressure,  given 
in  equation  (118),  reduces  to  the  following  second  order  (Fresnel) 
approximation  : 


P(r„) 


jk  P„ 

_ 0 

2tt  r 


rV- 


1 M 


dx  dv 


(121) 


jk  P e -^’^o 
■'  o 

2ir  r_ 


e-j*^  dx 

dy 

^1 

1 

>'l  1 

Substituting  u^  = , v^  = 

‘ o 0 
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the  integrals  in  Equation  (121)  become  the  Cornu  spirals  (see 
E.  Skudrryk , pp  629,  Equation  156): 


P (r^)  = j e'-^^^'o 

1_  ! 

p du 

1_ 

■''2  iv 

0 0 

n J 

/I  j 

''l 

where 


^1,2  = ^ 


1,2 


■ >'1.2 


Each  of  the  Cornu  spirals,  see  Figure  (38),  have  an  asymptotic 
value  of  (1-j)  if  the  limits  of  the  integral  are  equal  to  or 

exceed  ±— i—  . 

/I 


Thus,  in  the  limit  of  a large  sized  reflector,  the  reflected 
pressure  in  Equation  (1 22) approaches  the  value 


as  one  would  expect  for  reflection  of  a plane  wave  from  an  infinite 
plane  rigid  reflector 


If  the  field  point  projection  falls  on  the  center  of  the 
rectangular  reflector,  then  x,  , = ±a  and  y,  - = ±b . The 
limiting  case  of  large  plane  reflector  is  reached  when  the 
limits  u,  or  v,  exceed  — , i.e.  when 

n 


1_ 

/I 


and 


to 
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(2a) 
“7 — 


and 


(2b) 

X 


Equation  (123)  shows  that  the  linear  dimensions  of  the  reflector, 
normalized  to  the  acoustic  wavelength,  must  be  larger  than  the 
square  root  of  the  normalized  range  r^/X.  For  a square  planar, 

2a  = 2b  = h,  the  maximum  range  r^  where  the  geometrical  optic 
theory  applies  is  given  by 


2a 

T" 


2b 


h 

r 


m 


The  limit  range  r^^  is  shown  in  Figure  (41)  in  terms  of  the 
frequency  and  the  linear  dimension  of  a square  reflector  h. 

It  is  instructive  to  obtain  exact  values  of  the  integrals 
in  Equation  (122)  for  various  values  of  the  normalized  dimensions 
of  the  reflector.  If  2a  = 2b  = h,  and  if  the  field  point 
projection  is  located  at  the  center  of  the  square  reflector. 


then  tlie  amplitude  of  the  normalized  pressure 
given  by: 


P(ro)/Po 


IS 


P(ro) 


r\i 


e-j"'  dv 


u = 


The  exact  value  of  the  normalized  reflected  pressure  for  a square 

reflector  (h  x h)  is  shown  in  the  following  table  for  various 

locations  of  the  field  point  r„/r  : 

‘ o m 


(123) 


(124) 


(125a) 
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u 

r /r 
0 m 

Pfr  1/P  -Fresnel , 

0 Equation  (125a) 

PTr  1/P  .Fraunhofer, 

^ o-^'  0 Equation  (125b) 

0.10 

50.00 

0.02 

.02 

0.30 

5.50 

0.18 

0.18 

0.50 

2.00 

0.47 

0.50 

0.71 

1.00 

*1.00* 

0.50 

*1.00*  (Limit) 
1.22 

*1.00*  (Limit) 
2.00 

1.20 

0.35 

1.74 

3.30 

1.50 

0.20 

1.39 

5.00 

2.00 

0.10 

0.71 

10.00 

3.00 

0.05 

1.28 

20.00 

Thus  , in  the  range  r 


the  laws  of  geometrical  optics 


o m 

apply,  and  the  reflected  pressure  is  equal  to  the  incident 
pressure,  i.e.  the  reflected  pressure  is  the  same  as  if  the  wave 
was  geometrically  reflected  at  the  flat  surface  of  the  reflector, 
see  above  table.  The  condition  in  equation  (124)  is  exactly 
the  same  formula  as  the  one  given  in  equation  (113)  as  the 
requirement  for  the  existance  of  a stationary  phase  point,  i.e. 

7 

that  r < r = h /X. 
m 

For  distances  r > r , the  Fresnel  integral  should  reduce 
to  the  Fraunhofer  approximation.  For  a field  point  on  the 
normal  axis,  x,y  <<  r^  and  the  exponentials  can  be  replaced  by 
unity,  then  the  Fresnel  integral,  Fq . (122),  reduced  to  the 

expression  given  in  Equation  (120)  with  D (9,<4)  = 1: 


jk 


P^(4ab) 


Trf  r 


J 


Po(4ab) 


•jkr 


0 = j 


-jkr 
e •’  o 


m 


(125b) 


which  has  the  spherical  (1/r^)  falloff.  Equation  (125b)  predicts 
the  reflected  amplitude  accurately  up  to  , see  above  table. 
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This  result  is  equivalent  to  assuming  that  the  radiation 
is  geometrically  reflected  provided  the  radiation  makes  an 
angle  20  with  the  field  point  that  is  at  least  equal  to 
which,  for  a square  reflector,  see  Figure  (39), 


sin 


min 


min 


h/  2r 

m 


A/2h  , 


which  is  also  given  in  liquation  (114). 

If  the  incident  field  is  produced  by  a point  source  located 
at  field  point,  then  the  backscattered  field  can  be  computed 
by  the  use  of  Equation  (118).  Let  the  acoustic  field  of  the 
point  source  be  described  by 


where  r is  given  by  Equation  (117).  Using  Equation  (127)  in 
the  expression  for  the  reflected  pressure  in  Equation  (118), 
one  obtains  , with  rs=:  r^  , 


P(Tq  ,6  ,<f) 


^0 


■ + a 

fb 

exp  2jk 

^ -a  ■ 

-b 

cos  6 + y cos 


[x  cos  9 + v cos  q] 

27; 


dx  dy  . 


a b 

If  — <<  i and  — <<  1 , and  if  the  field-source  points  are  not 
o 

at  normal  incidence,  then  the  integral  can  be  simplified,  to  a 
first  order  (Fraunhofer)  approximation,  as  follows: 


(126) 


(127) 
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P(r,0,4i)  = 


jk  e 

2tt  r*" 
o 


ra 


-a  i 


p2jk[x  cos  0 + y cos  $] 


-b 


jk  p 


2tt  r‘ 


D(0,<.) 


(129) 


where  the  directiv'ity  factor  D is  giv'^cn  by 


t 


DC0,4))  = glg..-L“K^-.5:ps_Q) 

(2ka  cos  0) 


sin  (2kb  cos  4)) 
(2kb  cos  4>) 


sin  (ka  cos  0)  cos  (ka  cos  6) 
(ka  cos  0) 


sin  (kb  cos  f)  cos  ^kh  cos  :{)) 
(T^ cos  t J 


(130) 


Maximum  backscattering  is  attained  when  the  field-source  points  are 
at  normal  incidence,  i.e.  when  6 = <P  = j,  where  D — ► 1.0.  The 
backscattered  pressure,  at  large  distances  where  r >>  a,b, 
decreases  as  1/r^  . This  form  of  backscattering  is  equivalent 
to  spherical  spreading  from  source  to  reflector  and  again  from 
the  reflector  to  the  field  point.  The  phase  change  of  (2kr^) 
reinforces  this  fact.  For  such  large  distances,  the  laws  of 
geometrical  optics  do  not  apply. 

The  predeeding  computations  are  not  valid  for  distances 
closer  than  rj^ , as  defined  in  Equation  (111).  For  normal 
incidence,  0 = * tt/2  and  one  must  retain  more  terms  in  the 

eikonal  expression  in  Equation  (128).  Thus,  the  backscattered 
field  becomes  , 


jk  P 


‘2  -u  2 , 

/''o  dx 


ZTT  r 


^2  2 

-jky^/r  , 
e o dy 
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j P 
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du 


1_ 

/I 


r^- 


. TT  2 

•Jt  V 


dv 


where  u' 


2k  x' 
71  r. 


2k 
77  r 
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1,2 
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2 y 


1,2 


1 ,2 


/tt: 


The  Cornu  spiral  integrals,  see  Figure  (33),  have  an  asymptotic 

value  of  (1-j)  if  the  limits  are  equal  or  exceed  ± — 

/T 


In  the  limit  of  a large  rigid  reflector,  a ,b  >>  /X  r^  or 


“1,2 


PCr„) 


n -2jkr^ 
P e 0 
o 

2r_ 


viiiich  represents  the  reflected  field  of  a point  source  by  an  infinite  rigid 
reflector.  If  the  field-source  points  projection  on  the  plane  falls 
the  center  of  the  reflector,  then  Xj^  2 ~ ^1*^  limiting 

case  of  large  reflector  is  reached  when  u^^  t and  v^  ^ exceed 
± 1//T  , i . e . when 


2a 


> — and 


^0 

/I 

or 

(2a) 

/o 

V 7X 

and 

i|bi. 

2b 


/X  r. 


/2 


(131) 


(132) 

at 


(133) 
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Equation  (133)  shows  that  the  wave  length -norma  1 i zed  linear 
dimensions  of  the  reflector  must  be  larger  than  the  square  root 
of  half  of  the  wavelength -normalized  range.  For  a square  reflector, 
let  2a  = 2b  = h,  then  the  maximum  range  r , where  the  geometrical 
optic  theory  applies,  is  given  by 


2a 

T~  " 


(134) 


The  limit  range  r^  is  the  same  one  predicted  by  the  requirement 
for  the  existance  of  the  stationary  phase  point,  see  Equation  (111). 
If  the  field  point  is  not  on  the  axis  of  the  rigid  reflector,  we 
may  assume  the  reflector  to  consist  of  one  of  infinite  size  with 
two  (or  four)  pressure  release  reflectors  adjacent  to  it.  The 
resultant  effect  then  will  be  that  of  the  finite  reflector. 

The  above  computations  have  been  performed  as  if  the  reflector 
v/as  in  a baffle.  This  assumption  is  always  made  in  computing 
transducer  arrays.  It  leads  to  results  that  agree  well  with 
tests  on  arrays  for  all  angles  of  incidence,  except  those  close 
to  grazing. 

It  is  not  possible,  however,  to  deduce  the  results  for  a 
semi - inf inite  strip  from  the  proceeding  reflection  formulae 
because  the  strip  can  no  longer  be  approximated  by  a reflector 
in  a baffle.  For  a infinite  cylinder,  the  reflector  model 
will  represent  a reasonably  good  approximation.  The  exactly 
computed  reflection  curve  for  a infinite  cylinder  of  radius  R 
shows  that  for  a plane  wave  that  impinges  normally  to  its  axis, 
the  reflection  factor  of  a rigid  cylinder,  when  kR  ^ 1 , is 
approximately  equal  to  , (see  E.  Skudrzyk,  1971,  page  448). 

Thus  the  amplitude  that*^s  scattered  at  a rigid  infinite  cylinder 
is  : 


P(ro) 


(135) 
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A finite  cylinder  of  length  h one  would  expect  that 
the  cylindrical  spreading  (l//r^)  shown  in  equation  (135), 
would  be  valid  only  if  r^  < r^ , as  defined  in  equation  (124). 

In  the  far-field,  r^  > r^ , the  effect  of  the  finite  dimension 
of  the  cylinder  (h)  must  be  superimposed  on  the  curvature  effect 
given  in  equation  (135).  One  can  obtain  a good  approximation 
for  the  far-field  correction  by  multiplying  the  expression  in 
Equation  (135)  by  the  reflection  from  two  finite  edges  of 
length  h.  This  correction  factor  can  be  obtained  from  equation 
(122)  by  neglecting  the  second  integral,  resulting  in  a factor 


Thus,  the  total  reflected  field  from  an  insonified 


finite  cylinder  becomes: 


P(ro) 


Tr 


xr. 


= P 


h_ 

r 


r > r 
o m 


(136) 


where  r^  is  given  by  equation  (124).  expression  for  the 

reflected  field  from  a finite  cylinder  shows  the  spherical 
falloff  with  the  field  point  distance  r^ . 

I i,.. 

Similar  cone  lus ions -can  also  be  drawn  by  considering  the 

Fresnel  solution  for  a circular  reflector  or  by  analyzing  the 

exponent  in  the  diffraction  integral  near  a stationary  pliase 

point.  The  corresponding  integrand  always  contributes  fully 
2 > 

if  u = 1/2  (i.e.  if  r < r ):  where  u is  the  same  variable 
^ o m-^  ’ 

as  in  the  Fresnel  integrals  see  Equation  (122).  This  condition 
then  represents  the  limiting  range  of  the  reflection  or  diffraction 
field  that  is  satisfied  by  the  geometric  diffraction  theory 
in  which  the  field  amplitude  decreases  as  l/»'^  (see  proceeding 
section).  If  u"  < 1/2  (i.e.  if  r^  > r^)  , the  exponential  can 
be  replaced  by  1,  and  the  reflection  or  diffraction  field  decreases 
t/ith  the  distance  as  I/Tq.  as  if  the  reflector  was  a point  source. 
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Proper  accounting  must  be  made  for  the  distance  dependence 
when  bodies  of  finite  size  are  being  investigated.  Since  all 
finite  bodies,  in  the  far-field,  have  a 1/r^  distance  dependance  , 
the  one-dimensional  solution  (infinite  edge,  which  has  a l//r^ 

dependence)  when  applied  to  finite  edges  and  wedges,  should 
have  its  amplitude  multiplied  by  the  factor  / r^/ . Thus,  while 
the  amplitude  of  the  field  from  an  infinite  edge  falls  off  at 
a rate  of  3 dB  per  doubling  of  distance,  that  of  the  finite 
edge  falls  off  at  a rate  of  6 dB  per  doubling  of  distance. 

The  error  introduced  by  the  use  of  the  edge  diffraction  model, 
at  distance  r > r , compared  to  the  finite  diffractor  is 
illustrated  in  the  table  below  for  doubling  of  distance  relative 


<^^dge/^’finite)  dB  ^ 

1 

O 

o 

2 

3.0 

4 

6.0 

00 
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1 4 . DIFFRACTION  I-IELD  I-ROM  EDCIES  01-  FIN'ITi:  LEN'GTH 

Edge  diffraction  is  a highly  directive  phenomenon.  It 
has  been  shown  by  Fermat's  principle,  see  Appendix  HI.  that 
the  radiation  scattered  by  each  edge  element  ds  lies  on  a 
cone  whose  semi -apex  angle  is  the  same  as  the  angle  of  the 
incident  ray  with  the  element.  Thus,  the  scattering  phenomenon 
associated  with  edge  diffraction  is  somewhat  similar  to  that 
of  reflection  - the  angle  of  scattering  at  an  element  ds  is 
equal  to  the  angle  of  incidence.  However,  because  the  scattering 
element  is  not  a surface  but  a line  element,  scattering  occurs 
symmetrically  around  that  line  element,  i.e.  the  scattered 
rays  lie  on  the  surface  of  a cone.  For  example,  a p.  ae  wave 
that  is  incident  at  an  angle  S with  the  edge  has  a scattering 
cone  angle  6,  see  Figure  (42a).  For  a plane  wave  normally 
incident  on  an  edge,  i.e.  for  0 = tt/2  , the  scattered  rays  lie 
in  a plane  that  is  perpendicular  to  the  edge  element,  see  Fig.  fl2b) . 

For  a finite  edge  of  length  Z,  the  scattered  radiation 
is  confined  in  the  volume  between  two  cones  , whose  apexes  are 
at  each  end  of  edge,  S'^e  Figure  (42c).  Geometrically  for  a 
plane  wave  incident  at  an  angle  0,  the  scattered  field  can 
only  be  received  in  a region  corresponding  to  the  volume  of 
a conical  shell  of  thickness  equal  to  the  projection  of  the 
edge  length  on  wave  front.  If  the  receiver  is  in  this  conical 
shell  region,  then  the  apex  of  the  scattering  cone  originates 
from  a point  on  the  edge.  Mathematically,  this  point  corresponds 
to  the  stationary  phase  point  for  the  receiver  position.  The 
edge  must  be  long  enough  so  that  a stationary  phase  point  can 
be  formed  on  it.  If  there  is  no  stationary  phase  point  on  the 
finite  edge,  then  there  is  no  scattered  or  diffracted  field 
in  the  sense  of  the  geometrical  theory. 

Next  let  us  assume  a spherical  incident  wave  on  a finite 
edge  of  length  I,  see  Figure  (43).  The  condition  for  the 
stationary  phase  point  is  satisfied  by,  see  Figure  (17), 
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C137) 


Next,  let  us  examine  the  scattered  field  in  the  plane 
containing  the  source  L and  the  scattering  edge.  The  opening 
of  the  scattered  field  then  is 


H s in  B 

P 


(138) 


where  p is  the  distance  from  the  source  L to  the  projection 
point  P,  see  Figure  (43).  If  the  distance  p>>S.,  the  scattering 
opening  angle  e is  very  small . According  to  the  geometrical 
theory  of  diffraction,  the  scattered  field  at  receiver  points 
located  approximately  as  far  from  the  edge  as  the  source  would 
be  confined  to  a conical  shell  not  very  much  thicker  than  the 
twice  the  projection  of  the  edge  on  a wave  front.  This  means 
that  the  field  that  is  scattered  by  an  edge  of  a length  of 
15  feet  can  only  be  observed  over  a cross  section  of  not  much 
more  than  30  feet. 

This  means  that  the  geometrical  theory  of  diffraction 
cannot  be  applied  to  edges  of  finite  length.  In  this  cabO, 
the  geometrical  theory  gives  expressions  for  the  scattered 
sound  field  that  are  both  correct  dimensionally  and  in  absolute 
value,  provided  it  is  used  in  the  range  in  which  it  is  applicable, 

, Equations  (111)  or  (113).  In  the  far  field  [r  >r  , 
F.quations  (111)  or  (113)],  the  exact  diffraction  formula  for  a 
finite  sized  scatterer  are  very  similar  to  the  radiation 
formula  of  finite  baffled  array.  However,  an  even  more  accurate 
solution  for  the  diffraction  field  from  a finite  edge  can  be 
obtained  using  the  Rubinowicz-Kirchoff  integral  given  by  Equation  (21) 


d 


A 

W 


,-jk  (o  + r) 


cos  (n,r) 


sin  (p,ds)  ds  (139) 
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edge 


rp 


1 + (cos  r,p) 
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f 


where  the  factor  — ^ — represents  the  incident  spherical  wave, 

p and  r arc  the  distances  of  the  source  and  receiv^er  to  the  edge 
element  ds . 

Consider  a finite  curved  edge  of  length  I,  which  is  small 
compared  to  p and  r,  such  that  one  may  substitute  for  the  term 
1/ro  in  the  amplitude  factor  bv  l/r„p„  , where  r and  p are  the 
distances  r and  to  origin  s = 0,  located  at  the  center  of  the 
edge.  Retaining  only  r^  and  p^  in  the  amplitude  factor  l/rp, 
the  intergrai  , equation  (139),  simplifies  to: 


f c 
i edge 


jk(p  + r)  , cos  (n  , r)  sin  S ds 
1 + cos  (r,p) 


(140) 


where  £(sj  ■*  ) (r,ds)  is  the  angle  the  radius  vector  p of  the 
source  point  makes  with  the  edge  element  ds . For  a straight  edge, 
the  shadow  boundary  is  plane,  and  consequently  n = constant. 

The  influence  of  the  curvature  of  the  edge  is  contained  in  the 
functional  relations  of  r(s),  e(s),  B(s)  and  y(s),  where 
y(s)  = ) (r,ds)  is  the  angle  between  the  radius  vector  r of  the 

field  point  and  tlie  edge  element  ds  . The  sum  of  the  radii  r and 

p can  be  expanded  in  a Taylors'  series  on  the  distance  along 

the  curved  edge  s,  sec  Appendix  IV,  the  leading  term  of  the 

exponential  being  p^  + r^ . 

For  a fixed  position  of  the  source  and  field  points,  it 
is  convenient  to  establish  a cartensian  coordinate  system 
centered  at  the  diffracting  element  (ds) , see  Figure  (44)- 
Let  the  positive  z-direction  be  directed  in  the  direction  of 
the  incident  ray,  the  x-axis  in  the  plane  of  ds  and  u,  , (i.e.  in 
the  plane  of  the  shadow  boundary)  where  u ^ are  unit  vectors 
of  the  cartesian  coordinate  system.  The  direction  vectors  of 
the  field  and  source  points  are: 


P u. 


(141) 


-r 


X u^  - 


y Uy  ♦ 


z u. 


(142) 
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The  unit  tangent  vector  of  the  edge  is  given  by: 

= U - - (X  u,  ^ y Uy  ^ z (143) 

The  unit  normal  vector  of  the  shadow  boundary,  pointing  into  the 
shadow  zone,  is  then  given  by: 


n(s) 


t (s)  = y u^  - X u 


(144) 


The  various  angular  functions  of  the  integrand  in  Equation 
(140)  can  now  be  evaluated  in  terms  of  the  new  coordinate  system. 
Thus  , 


cos 

(n,r)  = 

-♦  r 
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(145) 
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(146) 

cos 

B = cos 

(0,ds) 

= cos 

(p,t) 

• 

“ ■“  z 

(147) 

cos 

Y = cos 

(r  ,ds) 

-*■ 

r 

t = 

XX  + yy  + zz 

(148) 

r 

r 

Combining  those  expressions  for  the  angular  function,  one 
obtains : 

cos  (n,r)  yx  - xy 

1 cos  (r,p)  ^ 


r 


^ V 

r ' 


(149) 


1 - - 


For  large  distance,  r>>l,  the  coordinates  of  the  numerator  can 
be  replaced  by  their  values  at  the  center  of  the  diffracting 
edge  (Xq  , y^ , z^) . Thus,  the  directivity  factor.  Equation  (149), 
reduces  to: 


cos  (n,r) 

1 ♦ cos  (r,o) 


r_ 


(150) 
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where  the  term  z/r  - has  been  neglected  when  compared 

to  the  other  terms.  It  must  be  noted,  however,  that  such  an 
approximation  to  the  denominator  of  Equation  fl49')  is  not  valid 
if  the  field  point  falls  on  the  shadow  boundary,  because 
cos  (r,p)  = -1.  Thus,  near  the  shadow  boundary  the  denominator 
of  Equation  (150)  must  be  left  as  it  is  given  in  Equation  (149). 

Introducing  a polar  coordinate  system  at  (ds)  , whose  z-axis 
coincides  with  z-axis  of  the  cartesian  system,  i.e. 


_z  = r cos  9 

y = r sin  6 sin  (151) 

X = r sin  9 cos  . 

The  angular  factor.  Equation  (150),  thereby  reduces  to: 

sin  9 sin  6 ^ - sin  9 cos  (b  v 
1 = 

1-  cos  9 

(152) 

cot  [x  sin  Cfi  - y cos  i] 


The  diffraction  integral  for  a finite  curved  edge  in 
Equation  (140)  becomes 

•112 


$ 


A 

^^0^0 


[x  sin  (}>  - y cos  iji]  ds 


e sin  g cot  fy 


(153) 


where  the  exponential  term  (p+r)  can  be  written  in  terms  of  the 
edge  variable  (s)  as  shown  in  Appendix  IV.  In  the  special 
case  of  a straight  edge  the  distance  from  the  source  and  field 
points  to  the  edge  may  be  approximated  by 

p * p^  - s cos  6 

r * r^  - s cos  6 
o 
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and  the  diffraction  integral,  Equations  (153),  can  be  written 
as 


$ 


d 


A e 


-jk(Do-^ro) 

^oPo 


rH/2 

J -Z/2 


^jksC''os  B + cos  y) 


sin  B cot  (j)  [x  sin  4 * y cos  4]  ds  (154) 

It  will  be  recalled  that  in  the  Kirchoff  theory  of  diffraction, 

the  orientation  of  the  plane  of  the  diffractor  is  unimportant  and 
that  it  is  only  the  discontinuity  introduced  by  the  edge  that  is 
of  interest.  Therefore,  if  the  edge  lies  in  the  y-z  plane  and  is 
inclined  by  an  angle  a = tt/2  - 6 to  the  x-axis;  we  obtain  the 

results  that  y = y = 0 and  x = cos  a = sin  B respectively.  In 

addition,  the  angular  functions  B(s)  and  y(s)  becomes: 


cos  B = 2 = constant 

X 2 

cos  Y " — sin  B - — cos  B 
r r 

o 0 

= sin  6 sin  B cos  4 - cos  9 cos  B (155) 

thus,  defining  the  sum  of  these  consines  as 

B = cos  B + cos  Y = sin  6 sin  B cos  6 + cos  B(l-cos  0)  (156) 


the  integral  for  a finite  straight  edge  in  Equation  (154)  can 
be  approximated  by: 

■ 1/2 
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0. 


sin^B  cot  (-2-)  sin  4 


.-"jk  B s jg 


■SL/2 


where  all  the  angles  were  assumed  to  be  constant  along  the  path  of 
integration,  since  l<<r,o.  The  final  result  of  the  integration  can 
be  shown  to  result: 


May  14,  1973 
EJS;SIH:ALS;  1j;is 


A Z sin  rkR£/: 

TFT  0 (TBITTn 


cot  (y)  sin^e  sin  o (158) 
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JtF  tf 


sin  CkB£/.l)  Q' 

n- R'.7 ’t 


(159) 


where 


0 • 2 „ 
— sin  3 


which  is  the  measure  of  the  acoustic  energy  incident  on  the  edge 
and  0'  = 9 + tt/2  is  the  polar  angle  as  measured  in  the  theory  of 
straight  edge.  The  factor  in  the  paranthesis  represents  the 

6 ** 

directivity  of  a dipole  line  array  of  length  Z,  and  the  tan 

term  is  the  directivity  factor  of  the  Kirchoff  straight  edge. 
Equation  (159)  shows  that  the  spreading  is  spherical  (the 
amplitude  being  proportional  to  Z./2r^)  and  that  the  directivity 
pattern  is  that  of  a linear  dipole  array  combined  with  that  of 
an  infinite  straight  edge.  Except  for  the  effect  of  the  angle 
functions,  the  diffracted  wave  is  the  same  as  if  the  incident 
wave  were  reflected  at  a sph._re  of  the  same  diameter  as  the 
length  of  the  diffracting  edge.  However,  the  last  formula  is 
valid  only  for  6>o  (i.e.  not  at  the  shadow  boundary)  because 
it  is  based  on  Equation  (150)  which  neglects  the  factor  (-'“o^/’^o 
in  the  denominator  of  the  integrand  as  given  by  Equation  (149). 

It  is  possible  to  derive  a solution  that  is  valid  very  near  the 
shadow  boundary  in  a similar  manner  as  this  has  been  done  for 
integral  in  Equation  (153).  However,  we  arc  not  really  interested 
in  this  solution  since'it  describes  basically  only  the  diffusion 
of  the  sound  field  across  the  the  boundary  of  the  geometric 
shadow.  Slightly  away  from  shadow  boundary,  Equation  (158)  is 
sufficiently  accurate. 

Equation  (158)  consists  of  two  groups  of  terms.  The  terms 
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not  in  the  brackets  are  the  familiar  diffraction  terms,  the 
cotangent  appearing  in  the  formula  because  9=0  counts  the 
derivation  in  angle  from  the  forward  direction.  The  bracketed 
factor  is  of  the  form  sin  x/x  and  is  the  same  as  the  directivity 
function  of  a linear  transducer  array.  For  field  points  not 
coincident  with  the  shadow  boundary,  this  factor  leads  to  maxima 
and  minima  that  are  similar  to  those  generated  by  a linear  array 
of  the  same  length.  The  amplitude  of  the  maxima  decreases 
proportionally  to 

I 

1 

2/2>  kr 

o 

where  B is  the  angle  factor.  Equation  (156),  that  on  the  average 
is  not  very  different  from  unity  and  is  the  amplitude  of 

the  diffracted  field  due  to  an  infinite  edge.  Thus,  Equation  (160) 
shows  that  for  a finite  edge  when  the  field  point  is  far  away, 
the  edge  is  not  long  enough  to  generate  a stationary  phase  point 
and  the  amplitude  decreases  from  the  diffracting  edge  at  a rate 
greater  (by  l//kr^)  than  for  an  infinitely  long  edge. 

It  should  be  noted  that  the  amplitude  of  the  maxima  further 
away  from  tbe  direction  of  the  main  beam  is  independent  of  the 
length  of  the  edge.  This  conclusion  may  seem  very  strange  at 
first,  but  a similar  relationship  can  be  obtained  by  considering 
linear  transducer  array  with,  for  instance,  a finite  number 
n of  elements.  The  magnitude  of  the  maxima  then  is  given  by 
(see  E.  Skudreyk,  1971,  page  598): 

^max  sin  A sin  [ (2v~  l)ir/2n)  ] (161 

where  is  the  directivity  function  for  the  maxima,  and 
A = (2v  l);r/2n,  v = 1,2  ...  n is  the  phase  difference 

between  successive  transducer  elements  in  the  direction  of  the 
maxima.  The  envelope  of  the  maxima  is  obtained  by  varying  v/n 
continuously  in  the  last  expression  from  0 to  n/Z.  The  amplitude 


n 


= A 


diff 


/Ttt  kr  B 
o 


(160) 
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of  the  main  maximum  is  therefore  proportional  to  the  number  of 
transducer  elements  in  the  linear  array. 

If  n is  great,  then  there  is  correspondingly  a great 
number  of  side  maxima  near  the  main  maximum  and  the  beam  is 
highly  directive.  As  a result,  little  of  the  array's  total 
energy  is  radiated  at  angles  away  from  the  main  beam.  If  the 
number  of  elements  is  small,  then  the  main  maximum  is  small 
too  and  much  of  the  energy  is  radiated  side-ways.  For  the 
same  transducer  element  spacing,  the  field  radiated  off  axis 
is  independent  of  the  number  of  elements  of  the  transducer 
array,  see  Figure  (45).  However,  there  is  one  essential  dif- 
ference; the  field  of  the  small  array  decreases  as  1/kr 
while  that  of  the  large  array  decreases  as  l//ITr  for  a 
much  greater  distance  from  the  array. 
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15.  SUMM\RY  AND  CONCLUSION'S 

The  acoustic  diffraction  and  hackscattering  fields  generated 
by  plane  diffractors  and  wedges  have  been  investigated  using  both 
the  geometrical  theory  of  diffraction  (GTD)  developed  by  Keller  and 
the  more  mathematically  formal  theory  based  on  the  K i rchof f -Rub inowicz 
integral  representation.  Comparisons  were  made  of  the  results  obtained 
using  the  two  methods.  A summary  of  those  results  and  their  physical 
interpretations  are  given  below'. 

15.1  The  Geometrical  Theory  of  Diffraction 

The  geometrical  theory  of  diffraction  (GTD)  is  based  on  the 
high  frequency  assumption  that  the  energy  of  a wave  can  be  confined 
to  tube- like  bundles,  and  propagated  along  ray  curves.  These  curves 
are  straight  lines  for  a homogenious  , isotropic  acoustic  media. 

Since  it  is  assumed  that  the  energy  is  equal  at  any  cross-section 
of  a bundle  of  the  rays,  the  geometrical  diffraction  field  is 
obtained  by  determining  the  variation  in  cross-sectional  area 
from  point  to  point  along  the  ray.  Curvature  effects  of  a 
diffracting  edge  are  shown  in  Section  (5)  to  be  obtained  directly 
by  this  method. 

Although,  the  GTD  gives  the  correct  distance  dependence 
(within  its  range  of  validity  for  the  size  of  the  diffractor), 
it  gives  no  indication  of  the  angular  variation  of  the  diffracted 
field.  Keller  determines  the  correct  directivity  factor  by 
forcing  his  solution  to  agree  with  the  exact  (Sommerfeld)  solution 
for  the  straight  edge  and  then  precedes  to  use  this  relation 
for  all  other  curved  edges.  (The  Ki rchof f- Rub inow icz  integral 
representation  directly  furnishes  the  proper  directivity  functions, 
or  at  least  a very  good  approximation  to  it.) 


15.2  The  Kircho f f - Rub inowic z Integral  Representation 


The  Ki rchof f - Rub inowicr  theory  of  diffraction  is  based  on 
the  solving  the  Helmhol  t z -iluyghen  integral  with  the  assumption 
that  the  diffractor  is  ideally  absorbent  and  therefore  no 
reflected  waves  occure.  Tliis  allows  the  diffracted  field  to 
be  investigated  independently.  (Reflections  may  be  accounted  for 
seperately  by  introducing  the  method  of  images,  see  Section  8). 

The  solution  of  the  resulting  Kirchoff  diffraction  field  is 
expressed  as  a surface  integral  over  the  shadow  boundary. 

Rubinowicz  reduces  this  surface  integration  into  a line  in- 
tegration about  the  perifery  edge  of  the  diffractor. 

The  resulting  edge  integral  can  then  be  evaluated  in 
various  kr  ranges  (k  acoustic  wave  number,  r is  the  distance 
from  the  edge).  It  is  interesting  to  note,  that  the  high  fre- 
quency limit  of  the  Kirchof f-R.ubinowicz  diffraction  field  is 
the  same  as  that  obtained  by  the  GTU ; i.e.  the  stationary  phase 
contribution  of  the  Kirchof f -Rubinowicz  integral  is  the  same 
as  the  direct  ray  contribution  of  the  GTD . 

The  Kirchof f-Rub inowicz  integral  representation  has  the 
advantage  of  being  mathematically  exact,  with  the  degree  of 
approximation  of  the  solution  (i.e.  the  diffracted  field)  always 
known.  This  method  provides  both  the  correct  distance  dependences 
and  angular  directivity,  and  is  therefore  more  useful  than  the 
GTD  for  curved  and  finite  edges. 

15.3  Diffraction  and  the  Gliadow  Boundary 

For  steady  state  phenomena  and  plane  diffractors,  it  generally 
appears  that  the  diffraction  field  can  be  attributed  to  the  discon- 
tinuity of  the  incident  wave  field  at  and  near  the  entire  shadow 
boundary.  The  shadow  boundary  acts  basically  like  a surface  layer 
of  dipole  sources.  The  effect  of  the  shadow  boundary  is  then  to 
generate  both  a wattless  near-field  and  a wave  which  propagates 
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from  the  region  near  the  edge  of  the  diffractor.  (This  is  similar 
to  the  radiation  field  of  a finite  plate  vibrating  below  coincidence; 
where  only  the  region  near  the  boundary  radiates  due  to  the 
uncancelled  volume  flow  between  adjacent  quarter  wavelength 
regions  of  the  plate.) 

The  near  field  contribution  arises  physically  from  the  dif- 
fusion of  momentum  from  the  incident  wave  into  the  shadow  region. 

This  near  field  is  correlated  with  the  incident  field  and  there- 
fore cannot  be  described  by  an  independent  energy  function.  In 
contrast,  the  field  generated  by  the  edge  area,  except  for  the 
region  in  the  immediate  vicinity  of  the  shadow  boundary,  is 
uncorrelated  with  the  incident  field  and  therefore  can  be 
represented  by  a seperated  energy  function  i.e.  the  diffracted 
field.  These  two  fields  establishes  tlic  continuous  transition 
of  the  pressure  field  from  the  insonified  to  the  shadow  regions. 

For  a transient  incident  acoustic  pulse,  such  as  a Dirac 
pulse,  the  orgination  of  the  diffracted  field  at  the  shadow 
boundary  can  be  demonstrated  directly. 

15.4  Diffraction  and  its  Flole  in  Determining  Backscatter ing 


Backscatter ing  from  a rigid  or  pressure  release  plane 
surface  is  made  up  primarily  of  two  phenoma : 

(1)  Spetral  reflections,  where  the  angle  of  the  incident 
wave  is  equal  but  opposite  to  the  reflected  wave,  and 

(2)  a diffraction  field  that  arises  at  the  edge  of  the 
surface . 

A concern  of  this  report  was  to  evaluate  what  effect  the 
diffraction  field  has  on  the  backscattered  field 

Two  diffraction  fields  arc  generated  as  a result  of  an 
edge;  one  is  caused  by  the  discontinuity  of  the  incident  field 
(i.c.  the  familiar  shadow  boundary),  the  second  by  the  dis- 
continuity in  the  reflected  field  due  to  the  termination  of  the 
reflecting  surface. 
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I The  diffraction  fields  in  both  cases  are  highly  directive 

f with  their  major  contribution  occuring  along  the  lines  of  the 

discontinuity,  in  the  direction  of  their  respective  propagation. 

Very  little  energy  is  propagated  in  the  opposite  directions. 

What  this  means  therefore,  is  that  the  major  contributions  to 
the  backscattcred  field  comes  from  tl\e  lobe  of  the  diffracted 
field  associated  with  the  reflected  wave.  Thus  the  application 
of  anechoic  coatings  to  a plane  reflector  will  reduce  the  amplitude 
of  the  reflected  wave  and  in  turn  tl\e  backscattering  caused  by 
the  lobes  of  the  diffracted  field  associated  with  the  reflected 
wave . 

Mathematically,  reflections  from  any  locally  reacting  surface 
may  be  included  in  the  Kirchoff  formulation  by  introducing  the 
method  of  images  as  is  explained  in  Section  C8) . 

15.5  Range  of  Validity  of  the  Geometrical  Theory  of  Diffraction 

The  GTD  has  been  shown  in  Section  (5)  to  be  based  on  Fermat's 
principle,  which  is  equivalent  to  the  necessity  of  having  a Kirchoff 
stationary  phase  point  located  well  within  the  length  of  the 
diffracting  edge.  The  diffracted  field  is  then  considered  to 
arise  from  the  edge  element  at  and  near  the  stationary  phase 
point.  For  edges  of  finite  length,  h,  the  range,  r^ , for  a 
stationary  phase  point  to  be  generated  (and  therclore  the  range 
of  validity  of  the  backscattering  obtained  by  the  GTUj  is 
given  by  r^  < N h /X  where  X is  the  acoustic  wave  lengtli  and 
N is  1 for  a plane  incident  wave  and  2 for  spherical  incident 
wave.  The  GTD,  therefore,  applies  only  within  a limited 
distance  range  from  the  diffractor.  Because  of  the  failure  of 
the  GTD  for  greater  distances  than  r^^^ , the  Kirchof  f - Rub  inowic: 
integral  representation  was  investigated  in  detail.  This 
representation  leads  to  solutions  of  the  diffracted  field  that 
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reliable  in  practically  the  whole  range  of  interest  and  lead 
to  the  correct  distance  dependence  regardless  of  both  the  size 
of  the  diffractor  and  the  frequency  (wave  1 engtli)  of  interest 
(Provided  that  the  diffractor  is  at  least  3 wavelengths  long). 

15.6  Diffraction  Field  from  a Finite  lidge 

Physically,  the  diffraction  field  from  a finite  edge  acts 
similar  to  the  radiated  field  from  a line  array  of  finite 
length.  Based  on  Ki rchof f -Rub inowic z theory  of  diffraction, 
the  diffracted  field  for  a finite  edge  was  shown.  Equation  (159), 
to  be  given  by  the  product  of  the  incident  acoustic  pressure, 
the  directivity  factor  of  a dipole  line  array,  and  the  directivity 
factor  of  the  Kirchoff  straight  edge.  This  relation  is  valid 
everywhere  except  near  the  shadow  boundary,  and  is  equivalent 
to  the  product  of  the  reflection  from  a sphere  of  the  same 
diameter  as  the  length  of  the  edge  and  some  combined  directivity 
factor . 
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APPENDIX  I 


Stationary  Phase  ^!ethod 

In  diffraction  work  it  is  frequently  necessary  to 
evaulate  integrals  of  the  form 

I #00 

I - S(s)  ds  = S(s)  ds 


where  c(s)  = (kr)  f(s),  k is  the  wave  number  and  r the 
range.  If  the  function  S(s)  is  a slowly  varying  function 
of  s in  the  range  of  integration  and  the  argument  of 
the  exponential  function  is  large,  i.e.  kr>>l,  the  integral 
may  be  evaluated  asymptotically  by  the  method  of  stationary 
phase.  Thus,  as  s varies  because  of  the  factor  kr,  e(s) 
will  vary  drastically  unless  it  passes  through  a maximum, 
a minimum  or  a point  of  inflexion  where  --j-g--—  = 0.  Since 

e(s)  represents  the  phase  of  the  integrand,  the  points  s^ 

where  --g|~ — vanishes  are  called  stationary  phase  points. 
Because  of  the  rapid  fluctuations  of  the  exponential 
factor  in  the  integrand  everywhere  except  near  stationary 
phase  points,  the  integrand  reduces  to  contributions  in 
the  neighborhoods  of  the  stationary  phase  points.  For 
example,  if  s^  = s^  is  a stationary  phase  point,  we  have 

8e(s  ) 

V _ A i 


Solving  this  equation  gives  s^,  and  if  there  are  more 
stationary  phase  points,  it  gives  the  roots  s^,  s^  ...  s^ 
for  all  the  stationary  phase  points.  The  exponential  e(s) 
may  then  be  expanded  into  a Taylor  series  near  this  stationary 
phase  point: 


e(s) 


3s 


Cs-s„) 


e"CSq)  + ... 
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3e(s  ) Cs-s  32^(5  ) 


(1-3) 


Higher  terms  need  not  be  considered  (unless  e"(sQ)  = 0) 
since  only  the  region  of  the  integration  in  the  vicinity 
of  the  stationary  phase  point,  (s-s^  — > 0),  contributes 
to  the  integral.  Substituting  Equation  (1-3)  into  (I-l) 
and  assuming  the  integrand  takes  on  the  value  at  the 
stationary  phase  point,  one  obtains  the  following 
asymptotic  value  for  the  integral: 


I - S(s^)  ej^^Cs^) 


gj(s-s^)2  c"(s^)/2  ds 


(1-4) 


Making  the  change  of  variable 


u = (s-Sjj)  /e''(s^)/2  , 


the  integral  may  then  be  written  in  the  form 
S(s  ) 

r\  ^ 


1 = 


+ "0 


~T 


du 


• i A 

The  integral  can  be  shown  to  equal  / it  e^^'  . Thus, 
value  of  the  complete  integral  is  thus  determined  by 
a single  value  of  the  functions  that  appear  in  the 
integrand,  and  is  given  by 


(1-5) 


I 


(1-6) 
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The  square  root  of  c"(s^)/2  has  been  taken  with  the  positive 
sign  because  we  have  assumed  that  ds  and  du  have  the  same 
sign  so  that  the  new  lower  and  upper  limits  of  integration 
is  again  and  respectively. 

At  a caustic  S(s^)  has  a pole.  The  stationary  phase 
method  then  must  be  generalized  also  for  this  integral 
by  separating  S(s^)  into  its  pole  contribution,  and  into 
its  regular  part.  (The  theory  for  this  situation  is  re- 
presented with  great  clarity  in  D.  S.  Jones,  1964,  page  689). 
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APPENDIX  II 

Derivation  of  the  Helmholts-Huyghens  Integral 

We  start  with  the  Green's  formula  relating  a volume  integral  to  a 
surface  integral 


- y 


dcr  = 


- W ^cp)  dx 


(Il-l) 


n = outside  normal  to  the  intergration  volume  x and  substitute  for  the 
field  'i  a solution  of  the  homogeneous  equation 


y^g  + k^g  = 0 


(H-2) 


that  is  finite  in  the  whole  space  of  integration.  The  right  hemd  side  of 
Eq.  (II-l)  then  is  zero.  To  obtain  the  field  at  a particular  point  r = 0, 
we  must  distinguish  this  point  from  all  other  points;  this  is  done  by 
identifying  g with  the  Green's  function  ” e for  a point  soiirce. 

The  region  near  r = 0 must  now  be  excluded  because  at  r = 0,  g(r)  is 
singular,  and  V^g(r)  is  undeterminate.  We  exclude  this  point  by  enclosing 
it  in  a little  sphere.  Hence  we  have 


0 = / da  (s\orface)  + / (little  sphere)  da 


^ 7 - f ^ i 

a r -*  0 

- • s - A I - ^ 

^ — A 


because 


and  r -» 0. 
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\ 


I 


Hence  we  have 


'P  •^)  da 
on 


To  obtain  this  integral,  the  field  point  r = 0 had  to  be  enclosed  by  the 
surface  of  integration.  If  the  field  point  is  outside  this  surface,  then 
the  last  integral  is  zero.  If  both  the  field  point  and  a source  are 
located  inside  the  surface  of  integration. 


P = P3_(Xty,z)  + 


(H-5) 


We  then  have  to  exclude  also  a small  sphere  around  = 0 from  the  space  of 
Integration,  and  obtain  also  a term  - A e 2 in  addition  to  the  term 
9 (O)  that  resulted  through  the  integration  around  r = 0.  The  sign  of 
this  source  term  is  negative  because  it  resxilts  from  the  first  term  of  the 
integrand  (which  has  a sign  that  is  opposite  to  that  of  the  second  term). 

We  thus  obtain,  as  would  be  expected 

q,(0)  + ^ J' (g^-P^)da  (H-6) 

0 

vhere  is  the  distance  of  the  source  from  the  field  point. 
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(II-7) 


The  above  resiilt  can  also  be  derived  purely  and  formally  with  the 
aid  of  the  spherical  5(r)  function  defined  by 


B(r)  4jtrdr  = Un  , ^ S(r)  Ujtrdr  = 0 

The  function  g then  satisfies  an  equation 


V^g  + = 6(r)  , (II-8) 

Recalling  that  cp  satisfies  the  homogeneous  equation  V^cp  + k^cp  = 0 , then 
on  the  right  hand  side  the  integral  becomes  (after  substituting  for 
V^g  ■ - k^g  + 5(r) ) 


^y^yyt<P(6(r)-  k^g)  - gV^cp]  dT  = J'fJ [<P6(r)  - g(^cp  + k^cp)]  dr  = 

ff r)dx  = tp(0) 


(II-9) 


o 


If  both  source  and  field  point  are  inside  surface  of  integration,  then 


+ k^g  = 6(r) 


V^tp  + k^cp  = BCrg) 


(II-IO) 


T 
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r 


^ g)<lT  = J" g[6(rg)  - k^cp]  - cp[6(r)  - k^g]  dr 

€ 

T 

= J'  [g8(rg)  - <p6(r)]  dT 

= ItJt  [g(rg)  - cp(0)] 


9(0)  = g(rg)  + ^ ff{g  da 

a 


t 


(Ii-u) 


(11-12) 
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APPENDIX  III 


Fermat's  Principle 


In  geometrical  optics  the  light  rays  are  defined  as  the  orthogonal 
trajectories  to  the  wave-fronts  cp(x,y,z)  = const,  by 


dr  , „ 

n — = grad  cp 
— as 


where  s is  the  distance  along  the  light  ray  and  n is  the  index  of 
refraction  of  the  medium  in  which  the  ray  travels. 

The  "ray  vector"  ns^  = ndr/ds  where  s^  is  its  unit  integral  vector 
and  satisfies  the  relation 


curl  ns*  = curl  (n  = curl  (grad  cp)  = 0 


and  by  integration  this  relation  over  an  area  a bounded  by  a closed  curve, 
we  obtain  the  Lagrange  relation; 


■/ 


curl  ns  • da  = P ns  • dr  = 0 
o / o 


This  means  that  the  Lagrange  integral, 


ns  • dr 
o 


is  Independent  of  the  path  of  integration.  By  comparing  the  Laigrange 
relation  over  an  element  of  the  actual  path  with  that  over  an  arbitrary 
curve  of  the  same  end  points  we  finally  arrive  at  Fermat's  principle 
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actual 

path 


nds 


nds 


arbitrary 

path 


(111-5) 


(Reference:  M.  Borne  and  E.  Wolf,  Principles  of  Optics,  Pergamon  Press, 

New  York,  I965). 

To  prove  Fermat's  principle,  we  take  a pencil  of  rays  and  compare  a 
segment  ^ with  an  arbitrary  curve  C joinging  P^^  and  P^ 

(Fig.  4 6).  Let  two  neighbor!  ig  orthogonal  trajectories  (wave-fronts)  of 
the  pencil  intersect  C in  and  and  C in  ^ and  Q^.  Further  let 
0^  be  the  point  of  intersection  of  the  trajectory  with  the  ray  C 

which  passes  through  Q^. 

Applying  Lagrange's  integral  relation  to  the  small  triangle  02^02^  ’ 
we  have 


{ris  . ‘ ~ ° 


(TIT-6) 


Now  from  the  definition  of  the  scalar  product 


(n?  . d?)^^  < . 


Further,  s is  orthogonsJ.  to  dr  on  the  wave-front,  so  that 


(ns*  . = 0 


Also,  since  corresponding  points  on  the  two  wave-fronts. 


i 
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On  substituting  from  the  last  three  relations  into  (lV-6) 


whence,  on  integration 


i 
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we  find  that 


(III-7) 


(III-8) 
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APPENDIX  IV 


Derivation  of  the  Rubinowicz  Integral 


The  surface  integral  Equation  (20)  in  the  text, 

can  bo  transformed  from  a surface  integral  over  the  face 
of  the  aperture  into  a line  integral  over  its  edge.  For 
this  transformation,  we  divide  the  surface  of  the  radiation 
cone,  shown  in  Figure  (47),  into  triangular  areas  with 
one  corner  at  the  apex  of  the  cone  and  two  sides  through 
the  end  points  of  the  edge  element  ds  of  the  aperture. 

The  element  ds  is  drawn  parallel  to  ds , and  has  a direction 
such  that  the  vector-cross  product 

da  = dp  X ds ' = dan 

represents  a unit  area  vector  of  the  shadow  boundary 
that  points  into  the  shadow  region.  Furthermore,  we 
have  (because  of  similar  triangles)  the  relationship: 

ds  * _ p 

a?-  ■ Ps  ’ 

and  therefore, 

da  = dp  X ds  ' = (dp  X ^ ds)  . 

^s 

To  evaluate  Equation  (20)  , we  still  need  the  distance 
r from  the  field  point  Q to  a point  on  the  shadow  boundary 
P.  This  distance  is  given  by  (see  Figure  47) 


I 

i 

I 


(IV-1) 


(IV-2) 


(IV-3) 


“ r * r^  (p  - Pg)  , 


(IV-4) 
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y > 


] 

-1 


f 


and  because  dp 


^ b g — ^ 

r • da  = r cos(n,r)da  = — — • (dp  x ds) 

= • (Pg  X ds)  dp  . 

Ps 

The  diffraction  integral.  Equation  (20)  of  the  text,  thus 
becoires 


d 


1 

J¥ 


(Pg  X ds) 


jk(p+r) 


where  the  integration  extends  from  p^  to  °°  , and  the  element 
ds  runs  over  the  edge  s of  the  aperture. 

The  p integration  can  be  performed  explicitly.  It 
is  expedient  to  introduce  the  vector  LP  = 2a: 

L?  = 2a  = LQ-PQ  = P-  r . 

If  we  square  this  equation  and  add  2r  • p to  both  sides, 
we  obtain 

rp  + r • p = 1/2  [ (r  + p)^  - 4a^]  ; 

which  when  differentiated  with  respect  to  dp  becomes: 
dr/dp  = cos  (r,p)  = r • p/rp 


(IV-5) 


(IV-6) 


(IV-7) 


(IV-8) 


(IV-9) 


t 
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Kith  the  aid  of  the  last  two  equations,  wo  easily  verify 
the  relation 

d fo  1 _ d [20  1 


rp  + r • pj 


(r  + p)  - 4a‘ 


,-jk(r  + p) 


ik  1 
I 1 • 


(IV-10) 


At  the  edge  of  the  aperture,  p = p^ , r = r^ , and  integration 
of  Equation  (IV-10)  leads  to  the  relation 


,-jk(r+p)  f-  ik  _ _1 


^ ^ do  - 


* “s’ 

'■s(''s»s*^s  • “s’ 


(lV-11) 


Hence  we  have  [see  Equations  (20)  and  (IV-6)]: 


♦ = a4> 


1 f ^ • r;  « d?] 

0 TrT  p r ^ 

Jedge  rp  > r . p 


g ■ W 


e ^ e ^ cos (r ,n) 
edge  ^ l + cos(r,p) 


sin  (p,ds)ds  , 


(IV-12) 


where  we  have  dropped  the  subscripts  s ; and  where  u = 1 
inside  the  cone  of  radiation,  while  a = 0 outside  the 
cone,  i.e.  in  the  geometric  shadow  region. 

The  integration  is  performed  around  the  contour  of 
the  aperture  in  a clockw'ise  sense,  when  looking  at  the 
light  source.  Note  that  n is  the  external  normal  to  the 
cone  through  the  shadow  boundary  (pointing  into  the  shadow 
region)  at  the  point  of  the  element  d s , and  Mr,p)  is  the 
angle  between  the  light  ray  from  the  source  and  the  radius 
vector  from  the  field  point  to  ds , as  shown  in  Figure  (7). 
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APPENDIX  V 


Frenet  Formulae 

Differentiation  of  position  vectors  locating  points 
on  curved  edges  are  known  as  Frenet  formulae.  Let  r 
be  the  position  vector  to  a curved  edge,  s be  the  coordinate 
system  for  the  edge  curve,  see  Figures  7 and  9. 

Letting  r = r(s),  r^  ^ r • r,  then 


where  t is  the  unit  tangent  vector  to  the  edge.  Likewise, 


(V-1) 


dr  /■-*•  r • t 

3s  = c°s(r,ds)  = 


(V-2) 


^ ^ (see  Figure  9) 

Ic 

where  ^ is  the  unit  vector  pointing  towards  the  center 
of  curvature  of  the  edge  at  s and  is  the  radius  of 
curvature  at  s.  It  can  therefore  be  seen  that 

idtt^-n 
t • t - 1 « • ati-  ‘ ° 

Thus,  differentiating  again  yields  the  relation 


(V-3) 


(V-4) 


7T  ^ Zs  r la?  j T [ a? 

ds  r 

. . ^ • i » c?  • S) 

» - i cos^  (r,ds)  i ^ cos  (r  ,ii) 

- i sin^(r,ds)  ^ cos  (r,H) 


(V-5) 
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Differentiating  the  unit  vector  h,  results  in 


dt  t 
" T 


(V-6) 


where  ? is  the  unit  vector  along  the  binormal  and  T is 
the  t rsion.  Thus, 


d'^r  3 ■^^2  dr  2 


C?  • 3?  - 3 


57  7 - r- 


dr 


2 I ^ 

3 tr  ■ t)  r • 35- 
r 


1 dr 

"7  3F 

r 


r“R 


1 

dr  t 

1 4.  1 

IT:  . d?;i 

■3 — • h 

[JF  J 

1 

Ji-J 

rR 


K) 


h) 


= ^ cos^  (r,ds)-  cos  cos  (r,ds) 


cos  (r,K)  - cos  (r,ds)  - ^ cos  (r,ds)  cos  (r  ,K) 


dR, 


+ cos  (r,^)  - cos  (r  ,h) 


(V-7) 


Rewriting  the  terms  in  CV-7),  one  can  show  that 


d^r 

57 


sin^  (r,ds)  cos  (r,ds)  - cos  Cr,ds) 


, . dR.  . ^ 

cos  (r,ii)  ♦ cos  (r,^)  - -j  -gj  cos  (r,h).  (V-8) 

k Rj^ 
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LEITNER 

WE«NER 

SPENCE 


«5  ^ 

0/2  T 


Figure  3-  Scattered  Pressure  Field 
on  the  Surface  of  a Rigid 
Disc,  Leitncr 

a.  Real  Part,  b.  Imaginary 
Part,  K refers  to  Kirchoff 
Solution,  c.  Absolute  Value 
by  Leitner,  approx.  by- 
Spence,  experimental  Data 
by  Weiner.  I and  S Refer 
to  the  Illuminated  and 
Shadow  Sides  of  the  Disc. 
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PLANE  OF 
DEFRACTING 
SCREEN 


DEFRACTING 
EDGE 


MAXIMUM  AMPLITUDE 


(1/2  THAT  OF 


UNDEFRACTED  FIELD) 


/V 


SHADOW 


SPACE 


SHADOW 

BOUNDARY 


Figure  8-  (a)  Section  through  diffraction  cone  containing 
incident  ray  and  shadow  boundary.  The  great 
intensity  of  the  diffracted  amplitude  near  the 
shadow  boundary  is  indicated  by  a thickening 
of  the  ellipse  on  the  drawing. 

(b)  Ar-litude  of  diffracted  rays  as  a function  of 
an^^'.e  from  shadow  region. 
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Figure  11-  A pair  of  neighboring  rays  incident  on  a curved 
edge  of  a flat  screen,  and  resulting  diffracted 
cones.  The  two  cones  of  diffracted  rays  intersect 
at  the  caustic,  which  is  at  the  distance  R from 
the  edge  along  one  of  the  diffracted  ravs.'^ 
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-i^ure  12-  The  focal  lines  of  a pencil.  A portion  of  the  axis 
of  rotation  and  a surface  of  revolution  whose  meri- 
dian section  is  the  cvolute  of  the  meridian  section 
of  the  wavefront.  CD.  S.  Jones,  The  Theory  of 
Electromagnetism,  Page  375). 


Figure  13-  Coordinates  for  the  Diffraction  from  the  Arc  of 
a Circular  Edge  of  Radius  a 
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Figure  15-  Shadow  Boundaries  of  Incident  and  Reflected 
(Image  Space)  Waves 


Figure  19-  Kirchoff  Solution  for  tlie  Diffraction  from  a Pressure  Release  Edge 
for  Various  Angles  of  Incidence  <|i 


1000.0 


Figure  20-  .Sommerf eld  (Exact)  Solution  for  the  Diff’-action  from  Rigid  Edge 


Sommerfeld  (Exact)  Solution  for  the  Diffraction  from  a Pressure  Release  Edge 
for  Various  Angles  of  Incidence  4> 


SOMMERFELD 

KIRCHOFF 


I60«,|200*  240®  280®  320®  360 


SOMMERFELD 

KIRCHOFF 


23-  Backs  cat  t ercd  I'ield  of  a Pressure  Release  lidge , 
Sominerfeld  and  Kirchoff  Solutions 


•I«I« 


:tion  from  a 30°  I’ressure  Release  Wedge,  for  Various  Angles  of  Incidence 


Figure  30-  Diffraction  from  a 90°  Pressure  Release  Wedge,  for  Various  Angles  of  Incidence 


Figure  31-  Diffraction  from 


091  = 


Figure  32-  Diffraction  from  a 120®  Pressure 


90°  and  120 
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(a) 


(b) 


Figure  39-  Limit  Range  and  Angles  of  the  Geometrical  Theory  of 
Diffraction  from  a Square  Scatterer  for  (a)  Point 
Source  and  (b)  Plane  have  Incidence 


LIMIT  OF  GEOMETRICAL  REFLECTION  AND  DIFFRACTION  ( r, 
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MAIN  MAXIMUM  LARGE 
NUMBER  OF  TRANSDUCERS 


MAGNITUDE  OFSIDE  LOBES  OF 
LARGE  ELEMENT  NUMBER  ARRAY 


1 1 

II!L 


SIDE  LOBES  OF  SMALL 
ELEMENT  NUMBER  ARRAY 


MAIN  MAXIMUM  SMALL 
NUMBER  OF  TRANSDUCERS 


Figure  45-  Radiation  from  a Finite  Linear  Array.  Ihe 

Off-Axis  Radiation  is  Independent  of  the  .Number 
of  Elements  in  the  Array. 
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Copy  No . 3 

Commander 

Naval  Ordnance  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  20360 
Attention:  Edward  Liszka,  Code  034A1 

Copy  No.  4 

Commander 

Naval  Ordnance  Systems  Command 
Department  of  the  Na^'y 
Washington,  D.  C.  20360 
Attention:  Edward  McKinney,  Code  034B 

Copy  No.  5 

Commander 

Naval  Ordnance  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  20360 
Attention:  William  August,  Code  0541A 

Copy  No . 6 

Commander 

Naval  Undersea  Center 
Pasadena  Laboratory 
3202  E.  Foothill  Blvd. 

Pasadena,  California  91107 
Attention:  Douglas  Chabries 

Copy  No . 7 

Commander 

Naval  Undersea  Center 
Pasadena  Laboratory 
3202  E.  Foothill  Blvd. 

Pasadena,  California  91107 
Attention:  Carl  Johansen 

Copy  No . 8 
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Commander 

Naval  Undersea  Center 
Pasadena  Laboratory 
3202  L.  Foothill  Blvd. 
Pasadena,  Calif.  91107 
Attention:  Joel  Young 

Copy  No . 9 


Commander 

Naval  Undersea  Center 
Pasadena  Laboratory 
3202  E.  Foothill  Blvd. 

Pasadena,  Calif.  91107 
Attention:  Jon  Reeves 

Copy  No.  10 

Commander 

Naval  Ordnance  Laboratory 
White  Oak 

Silver  Spring,  Maryland  20910 
Attention:  Guillermo  C.  Gaunaurd 

Copy  No.  11 


Commander 
Naval  Ordnance 
White  Oak 
Silver  Spring, 
Attention: 

Copy  No.  12 


Laboratory 

Maryland  20910 
Iter  Madigosky 


Commander 

Naval  Underwater  Systems  Center 
Newport  Laboratory 
Newport,  Rhode  Island  02840 
Attention:  Fred  Cancilliere 

Copy  No.  13 


Commander 

Naval  Underwater  Systems  Center 
Newport  Laboratory 
Newport,  Rhode  Island  02840 
Attention:  Dwain  Hartge 

Copy  No.  14 

Commander 

Naval  Underwater  Systems  Center 
Newport  Laboratory,  Bldg.  679 
Newport,  Rhode  Island  02840 
Attention:  Edward  Barile 

Copy  No.  15 


Continued 
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Commander 

Naval  Underwater  Systems  Center 
New  London  Laboratory 
New  London,  Connecticut  06321 
Attention:  Jan  Holland 

Copy  No.  16 

Commander 

Naval  Coastal  Systems  Laboratory 
Panama  City,  Florida  32401 
Attention:  Bruce  Nolte,  Code  762 

Copy  No.  17 

Commander 

Naval  Ship  Research  § Development  Center 
Washington,  D.  C. 

Attention:  Maurice  Sevik,  Code  900 

Copy  No.  18 

Defense  Documentation  Center 
Camron  Station 
Alexander,  Virginia  22314 
Copy  No.  19-30 

Via:  Commander  (ORD-632) 

Naval  Ordnance  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C.  20360 


